An-All v  VU6  DELAWARE  UNI*  NEWARK  OCPT  OP  CIVIL  ENGINEERING  F/fl  11/3 

■AWE  AINCRATICS  ANO  SEDIMENT  SUSPENSION  at  VAVC  0REAKIN8  POINT. (U) 
JUN  82  P  A  HU  ANO  *  H  NANO  N000 1  **-8 1  -K-0297 

U*ICLASSIFI£0  CE-82-27  NL 


WAVE  KINEMATICS  AND  SEDIMENT  SUSPENSION 


AT  WAVE  BREAKING  POINT 


By 

Paul  Anche  Hwang  and  Hsiang  Wang 

Technical  Report  No.  13 
contract  No.  N00014-81-K-0297 
with  the  Office  of  Naval  Research,  Geography  Programs 


RESEARCH  REPORT  CE-82-27 

June  1982 


Ocean  Engineering  Program 


Department  of  Civil  Engineering 
University  of  Delaware 
Newark,  Delaware 


19711 


ABSTRACT 

4/ 

The  flow  field  and  sediment  suspension  at  wave 
breaking  piont  were  studied  in  the  laboratory. 

The  flow  field  concerns  mainly  the  drift  velocity 

on  sloping  beach.  The  drift  velocity  profiles  as  well  as 

! 

the  mean  drift  velocity  strength  were  established  through 
laboratory  measurements.  Various  existing  theories  were 

examined  and  compared  with  laboratory  results. 

/  -  . 

A  sediment  suspension  model  was  proposed  based  on 
the  diffusion  model  including  -thereon side rat  ions  of  the 
fall  velocity  reduction  in  an  oscillatory  flow  field,  the 
sediment  grading^and  the  degree  of  agitation  in  the  flow 
field.  The  model  reduced  the  dependence  of  the  data— 
fitted  coefficients  to  a  minimum.  Comparison  of  this 
model  with  data  from  various  sources  was  made.  The 
results  are  encouraging. 


lv 


The  breaking  wave  spectra  in  a  finite  water  depth 
were  discussed.  A  modified  Wallops  spectral  model  using 
finite  depth  Stokes  waves  and  cnoidal  or  solitary  waves 
for  different  ranges  of  water  depths  was  suggested.  The 
model  was  also  extended  by  using  the  stream  function  wave 
theory  which  was  found  to  give  the  most  satisfactory 
results  on  the  wave  spectral  prediction  both  in  terms  of 
the  variance  and  the  spectral  shape. 
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CHAPTER  1  INTRODUCTION 


It  is  a  common  practice  to  treat  the  breaking 
point  as  a  natural  boundary  for  various  studies  of  surf 
zone  problems.  The  breaking  point  itself,  however,  has 
not  been  widely  investigated.  In  an  attempt  to  obtain  a 
clearer  understanding  of  this  region,  the  author  studied 
the  Lagrangian  drift  velocity,  the  vertical  and 
horizontal  velocities,  the  wave  forms  and  the  sediment 
suspension  at  this  specific  location.  Due  to  the  diverse 
nature  of  the  problems  involved,  the  literature  review 
and  theoretical  derivations  of  each  topic  were  included 
in  respective  chapter  that  deals  with  the  specific  topic. 

Two  sets  of  laboratory  experiments  were  conducted 
for  the  research.  Chapter  2  deals  with  the  drift 
experiment  and  its  results.  Chapter  3  descibes  the 
sediment  experiments  which  collected  simultaneous  data  of 
surface  elevation,  vertical  and  horizontal  velocities  and 


Island  of  Sylt,  West  Germany  are  also  discussed  in  that 
chapter.  A  sediment  suspension  model  based  upon 
turbulent  diffusion  is  presented  for  the  core  region  of 
the  wave  field.  The  comparison  of  the  model  results  with 
present  experiments  as  well  as  experiments  by  other 
investigators  were  discussed. 

Chapter  4  studied  the  wave  and  velocity  spectra. 
A  modification  of  Wallops  spectral  model  using  finite 
depth  Stokes  wave  theory,  cnoidal  theory,  solitary  theory 
and  stream  function  wave  theory  is  proposed  for  finite 
depth  water  application.  However,  the  direct  spectra  of 
the  stream  function  generated  waves  are  found  most 
satisfactory  both  in  terms  of  the  total  variance  and  the 
spectral  shape.  The  comparison  of  wave  and  velocity 
spectra  of  various  wave  theories  with  those  of  the 
collected  data  are  presented. 

The  last  chapter  presented  the  conclusion  for  the 
entire  study  and  recommendations  for  future  studies. 


)  CHAPTER  2  MEAN  DRIFT  VELOCITY 


2.1  Introduction 

Drift  velocity  in  a  wave  field  is  the  velocity  of 
a  mean  current  resulting  from  open  orbital  motions  of 
water  particles.  It  plays  an  important  role  in  material 
transport  in  the  nearshore  zone  and  is  also  known  as  the 
mass  transport  velocity.  Stokes  (1847)  was  the  first  to 
point  out  this  property  and  treated  the  problem 
analytically  for  a  channel  of  infinite  length  and 
constant  depth.  His  solution,  based  upon  classic 
second-order  wave  theory,  resulted  in  a  mean  flow  in  the 
direction  of  wave  propagation. 

Longuet-Higgins  (1953)  treated  the  problem  in  a 
finite-length  channel  with  due  consideration  of  fluid 
viscosity  in  the  surface  and  bottom  boundaries.  His 
solution  indicated  that  the  drift  velocity  is  in  the 
direction  of  wave,  propagation  near  the  surface  and  the 
bottom  but  against  the  wave  in  the  middle  section.  Later 
on,  Huang  (1970),  Wang  and  Liang  (1975),  Mei  et  al. 
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(1972)  and  Dalrymple  (1976)  carried  out  further 
analytical  studies,  all  of  them  dealing  with  finite 
channel  of  horizontal  bottoms.  The  drift  velocity 
distributions  obtained  by  them,  except  Dalrymple' s  are 
all  similar  to  that  of  Longuet-Higgins.  Dalrymple' s 
formulation  was  based  on  Dean’s  (1965)  stream  function 
theory  and  his  solution  was  in  the  Eulerian  sense  as 
opposed  to  the  Lagrangian  drift  velocity  used  by  the 
other  investigators.  The  Lagrangian  drift  velocity  using 
stream  function  theory  is  developed  in  this  study.  The 
drift  velocity  profile  is  similar  to  that  of  the  Stokes. 

Experimentally,  Russell  and  Osorio  (1957)  and  Mei 
et  al.  (1972)  measured  drift  velocity  distribution  in  a 
horizontal  finite-length  channel.  Bijker  et  al.  (1974) 
investigated  the  influence  of  bottom  slope  on  drift 
velocity.  Most  of  their  measurements  were  made  before 
wave  breaking  with  emphasis  on  near-bottom  flow.  For  the 
bottom  drift  velocity,  Bijker  et  al.  modified 
Longuet-Higgins'  solution  by  the  incorporation  of 
shoaling  effect.  The  theoretical  values,  however,  are 
found  to  be  considerably  larger  than  those  obtained  in 

I- 

the  experiments. 
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The  present  wave  tank  study  concentrated  on 
measurements  of  the  drift  velocity  at  the  breaking  point 
under  different  types  of  breaking  waves  on  a  rigid  plane 
beach.  The  difference  in  breaker  types  would  probably 
influence  the  vertical  drift  velocity  profile  at  the 
breaking  point,  which  would  in  turn  control  the  sediment 
transport  pattern  across  the  surf  zone.  An  attempt  is 
made  here  to  clarify  the  drift  veloctity  patterns  under 
these  breakers.  The  data  are  compared  with  existing 
theories . 


Further  comparison  with  theories  is  made  for  a 
depth  averaged  offshore  drift  velocity  (defined  in  Eq. 
2-3).  The  smoothing  effect  of  the  vertical  integration 
eliminates  some  'of  the  data  scatter  due  tc  measurement. 
The  theories  thus  agree  far  better. 

For  convenience  of  reference,  we  define  (1)  the 
breaking  point  as  a  starting  point  of  wave  breaking, 
which  is  characteriszed  by  the  existence  of  maximum 
height  of  waves,  usually  showing  the  initiation  of  bubble 
and  foam  formation,  and  (2)  the  plunging  point  as  a  point 
where  the  shape  of  breaking  waves  completely  disintegrate 
when  their  crest  impinges  against  water  (see  Pig.  2.1). 


The  drift  velocity  distributions  were  largely  measured 


ging  breakers 
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within  these  limits.  The  material  of  this  chapter 
appeared  -partially  in  Sunamura',  et  al.  (1980)  and  Wang 
et  al.  (1982). 

2.2  Laboratory  Experiment 

2.2.1  Experimentation 

The  experiment  was  conducted  in  a  steel  wave 
flume  26  m  long,  0.5  m  wide  and  1.3  m  deep.  At  one  end 
is  a  Scotch  Yoke  piston  type  wave  generator.  A  glass 
section  7  m  long  is  located  10  m  from  the  generator.  A 
plane  beach  made  of  plywood  coated  with  black  paint  was 
placed  at  the  other  end  of  the  tank.  Two  beach  slopes, 
1/10  and  1/15,  were  used  during  the  test.  Sixteen  runs 
of  experiments  were  conducted  with  conditions  listed  in 
Table  2.1.  The  last  three  cases  were  difficult  to  obtain 
velocity  profiles  for,  but  the  wave  data  are  used  in 
Chapter  4  for  wave  form  analysis.  For  the  rest  of  the 
chapter,  the  discussion  is  concerned  with  the  first 


thirteen  runs. 

Based 

upon 

the 

relative  water  depth 

criterion,  i.e. 

,  the 

ratio 

of 

still  water  depth  to 

wavelength  at  breaking,  (h/L)^ ,  the  test  conditions  were 
all  in  the  lower  end  of  the  intermediate  depth  range 
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Table  2.1  Laboratory  test  conditions  of  drift  experiment 


case 

T 

Hc 

hc 

Hb 

U 

s 

kh 

(sec) 

(cm) 

(cm) 

(cm) 

(cm) 

(cm/s ) 

type* 

1 

1.00 

8.0 

40.0 

9.5 

12.1 

1.0 

1/15 

sp  .760 

2 

1.33 

11.7 

40.0 

14.0 

18.3 

2.2 

1/15 

sp  .694 

3 

1.34 

9.0 

40.0 

11.5 

13.1 

2.0 

1/15 

tr  .570 

4 

1.47 

11.9 

40.0 

14.5 

18.6 

2.3 

1/15 

tr  .625 

5 

1.56 

8.4 

40.0 

12.5 

13.7 

2.6 

1/15 

pi  .495 

6 

1.65 

13.6 

40.0 

16.5 

18.3 

2.6 

1/15 

pi  .545 

7 

1.89 

12.5 

40.0 

16.5 

17.1 

4.8 

1/15 

pi  .454 

8 

1.22 

10.0 

35.0 

11.7 

14.0 

0.9 

1/10 

pi  .657 

9 

1.50 

11.5 

45.0 

15.4 

16.7 

2.0 

1/10 

pi  .576 

10 

1.58 

11.0 

35.0 

12.4 

13.0 

2.4 

1/10 

pi  .475 

11 

1.61 

9.5 

35.0 

12.6 

15.0 

0.9 

1/10 

pi  .502 

12 

1.80 

8.0 

35.0 

12.4 

13.5 

.  1.5 

1/10 

pi  .422 

13 

1.89 

10.0 

35.0 

15.2 

15.1 

1.6 

1/10 

pi  .425 

14** 

0.84 

9.3 

40.0 

7.5 

14.0 

1/10 

sp  .972 

15** 

1.16 

12.2 

40.0 

13.0 

19.7 

1/10 

Sp  .818 

16** 

1.81 

6.4 

40.0 

9.5 

13.2 

1/10 

pi  .407 

17+ 

1.14 

11.6++ 

+ 

11.9 

15.3 

1/35 

sp  .749 

18+ 

2.21 

11.5++ 

+ 

11.7 

14.9 

1/35 

pi  .379 

*  sp:  spilling  **  no  mean  velocity  profiles 

tr:  transitional  ++  deep  water  wave  height 
pi:  plunging  based  on  calculation 

+  from-  Flick  et  al.  (1981),  no  mean  velocity  profile, 
constant  depth  between  45  to  55  cm 
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close  to  shallow  water  condition  (Fig.  2.2).  As  far  as 
breaker  types  are  concerned,  -  twd  of  the  tests  were 
spilling,  nine  of  them  were  plunging  breakers  and  the 
remaining  two  were  transition  types.  These  groupings 
were  based  upon  the  magnitudes  of  the  surf  zone  parameter 
as  suggested  by  Galvin  (1968),  Battjes  (1974)  and  later 


amplified  by 

Wang 

and  Yang  (1980).  The 

surf 

zone 

parameter  is 

a  non-dimensional 

quantity 

defined 

as 

H  /g  Ts ,  where 

H  is 

the  breaking 

wave  height, 

T 

is 

the 

wave  period, 

s  is 

the  bottom 

slope  and 

g 

is 

the 

gravitational  acceleration.  The  values  of  this  parameter 
for  the  test  conditions  and  the  breaker  type  designations 
are  shown  in  Fig.  2.3. 

Input  wave  characteristics  were  recorded  by  a 
capacitance  type  wave  gauge  installed  at  a  place  of 
constant  water  depth  1.5  m  in  front  of  the  beach  toe. 
The  wave-induced  water  particle  movement  near  the 
breaking  point  was  tracked  by  neutrally  buoyant, 
polystyrene  beads.  these  were  yellow-colored  with  a 
diameter  of  approximately  2  mm.  Not  all  the  beads  had 
exactly  the  same  density  as  the  water  which  actually 
resulted  in  a  more  uniform  vertical  distribution  of  beads 
in  the  water  column. 
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Test  range  in  terms  of  relative  water  depth  criterion 


PP  |  s=1/10 
i 

S=1/15  fji  (JT  s  s 
FLUNGING  1  SPILLING 


Breaker  types  in  terms  of  surf  parameter. 


Fig.  2.4  illustrated  the  apparatus  at  the 
measuring  -  section.  The  neutrally  buoyant  particles  were 
mixed  with  water  in  a  tracer  feeder  made  of  a  transparent 
plastic  tube  with  an  outer  diameter  of  1  cm.  One  tapered 
end  of  the  feeder  was  placed  on  the  bottom  under  the 
breaking  point,  while  the  other  end  was  connected  to  a 
rubber  ball  by  a  rubber  tube.  By  squeezing  the  ball 
during  wave  action,  many  tracers  were  ejected  from  the 
feeder  tip  into  the  water.  Some  tracers,  slightly 
lighter  than  water,  went  up  slowly,  so  that  almost 
uniform  tracer-distribution  around  the  breaking  point  was 
attained  in  a  few  wave  cycles  after  the  ejection. 

The  movement  of  the  tracers  was  filmed,  over  10 
or  more  consecutive  wave  cycles  for  each  run,  by  using  a 
16  mm  movie  camera  (Bolex  H-16  EBM  Electric)  operated  at 
a  speed  of  24  and  25  frames/sec  for  1/15  and  1/10  beach 
slopes,  respectively.  Color  films  (Kodachrome  40)  were 
used.  A  100  mm  telephoto  lens  was  applied  to  reduce  the 
parallax. 

The  trajectories  of  beads  were  traced,  as  many  as 
practical,  from  the  films  with  the  aid  of  a  film  analyzer 
(L-W  International,  Model  224A).  Examples  of  trajectory 
diagrams  are  shown  in  Fig.  2.5.  Greater  turbulence 
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Schematic  diagram  showing  measuring  apparatus. 
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caused  by  water-mass  impinging  at  the  plunging  point, 
particularly  in  the  case  of  plunging  breakers  of  larger 
heights,  affected  the  two-dimensional  icy  of  the  flow 
field  in  the  vicinity  of  the  breaking  point,  thus 
influencing  the  movement  of  the  tracers.  Consequently, 
long  term  trajectory  tracing  was  difficult  in  the 
plunging  breaker  cases  (compare  Figs.  2.5a  and  2.5b). 

Fig.  2.6  shows  the  definition  of  mean  drift 
velocity,  U,  used  in  our  data  analysis  at  a  point,  (x/2, 
y/2):  U=(a+b)/2T,  where  T  is  the  wave  period,  a  and  b 
are  the  displacements  of  a  tracer  measured  with  positive 
values  in  the  offshore  direction,  and  x  and  y  are 
horizontal  and  vertical  dimensions  of  the  trajectory 
under  consideration,  respectively.  Since  the  mean  drift 
velocity  obtained  from  the  trajectory  diagrams  showed 
considerable  spatial  scatter  (Fig.  2.7  as  an  example, 
the  variance  of  the  scatter  is  about  40%),  local 
averaging  was  done  by  grouping  the  data  as  demonstrated 
by  dashed  closures  in  Fig.  2.7.  Some  ambiguity  remained 
in  the  selection  of  averaging  areas. 


2.2.2  Laboratory  Results 


Tracer 

Trajectory 


Definition  of  mean  drift  velocity,  U 


On  the  basis  of  the  locally  averaged  values, 
isolines  were  drawn  to  indicate  drift  velocity 
distribution  near  the  breaking  point.  The  result  is 
illustrated  in  Fig.  2.6,  in  which  offshore  drift  is 
chosen  as  positive. 

The  lack  of  data  points  near  the  bottom  is 
partially  due  to  the  difficulty  of  differentiating  the 
tracers  crowded  in  the  vicinity  of  the  bottom.  The  film 
analyses  showed  that  there  always  exists  a  thin  bottom 
layer  having  onshore  drift  velocity,  irrespective  of 
breaker  types.  This  layer  is  about  1  cm  or  less  in  depth 
for  the  cases  of  1/15  slope,  and  may  exceed  3  cm  in  1/10 
beach. 

A  wide  mid-layer  shows  offshore  mass  transport 
(Pig.  2.8).  Since  continuity  has  to  be  satisfied, 
onshore  mass  transport  must  exist  near  the  water  surface 
to  balance  the  net  offshore  flow  in  the  water  column, 
i.e.,  the  mid-layer  offshore  discharge  minus  the  bottom 
layer  onshore  discharge.  Adequate  measurement  of  the 
upper  layer  onshore  mass  movement  velocity  was  not 
practical  by  the  tracer  method  because  a  large  quantity 
of  air  bubbles  formed  at  the  initial  stage  of  wave 
breaking. 


8. 
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Fig.  2.9  shows  some  of  the  vertical 
distributions  of  drift  velocity  at  the  breaking  point. 
The  values  were  based  on  isolines  intersecting  with  a 
vertical  plane  at  the  breaking  point.  Again, 
quantitative  determination  of  the  velocities  on  the 
bottom  and  at  the  surface  could  not  be  mae. 

2.2.3  Vertical  Distributions 

Vertical  distribution  of  drift  velocity  at  the 
breaking  point  is  qualitatively  similar  to  that  in  the 
offshore  zone  on  a  sloping  beach  (Bijker  et  al.,  1974) 
and,  to  an  extent,  similar  to  that  on  a  horizontal  bottom 
(Russel  and  Osorio,  1957;  Mei  et  al.,  1972).  Fig.  2.10 
compares  the  shapes  of  the  vertical  ditribution  as 
obtained  by  experiments  with  the  theories.  The  offshore 
drift  velocity  in  the  interior  region,  in  general,  shows 
a  fairly  uniform  vertical  distribution  much  more  uniform 
than  in  the  offshore  region.  These  smoother  profiles 
near  the  breaking  point  may  be  attributed  to  stong 
turbulent  mixing.  There  seems  to  be  little  difference  in 
the  offshore  drift  velocity  profiles  between  plunging  and 
spilling  breakers.  The  former  might  have  a  slight  bulge 
in  the  lower  or  middle  part  of  the  profie,  while  the 


Fig.  2.10  Comparison  of  shapes  of  drift  velocity  profiles 
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later  has  none. 

2.3  Mean  Offshore  Drift  Velocity 

At  present,  there  are  a  number  of  theories  which 
can  be  used  to  estimate  the  mean  drift  velocity.  All  cf 
these  are,  however,  for  water  of  constant  depth.  The 
results  are  summarized  here  with  reference  to  the 
coordinate  system  shown  in  Fig.  2.11. 

Lonauet-Hiqqins'  drift  veloc i ty 

Lonauet-Kiggins '  solution  (1953)  is  for  the 
interior  region  in  a  channel  of  constant  depth  with  a 
laminar  boundary  layer: 

U(z)=  f  Ufa,z/h)  (2-1) 

-tSinh 


where 


f ( kh,z/h)*{2cosh  2kh( z/h+1 ) +3+kh[ 3 (z/h)2 +4(z/h)+l] 
sinh2kh  +  3[2-5^~^  +  3/2H  (z/h)2  -1]}  (2-2) 


and  U(z)  is  the  drift  velocity  at  depth  z. 


We  define  here  the  depth  averaged  offshore  drift 
velocity  as 

(hi 

U=-l/h  ,  U(z)dz  (2-3) 

Ai 

where  h,'  and  h^  are  the  upper  ana  lower  limits  of  the 
offshore  flow,  respectively.  The  non-dimensional  drift 
velocity  becomes 

U/aV  k=F(kh)/4sinh2kh  (2-4) 

(hi. 

where  F(kh )*~Jly  f  ( kh,z/h)d(z/h) ,  h,  and  h*  are  the 
dimensionless  limits. 

For  large  values  of  kh  (deep  water  case),  h,=-l, 
h2»-l/3,  and  we  have 

U/aa<rk=2kh/27*0 . 074  kh  (2-5) 

for  shallow  approximation,  it  becomes 

U//gTi*0.060  k.4  (2-6) 


where  K*(H/h)t  . 
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Solution  by  Wang  and  Liang 

Wang  and  Liang  (1975)  derived  a  mass  transport 
solution  similar  to  that  of  Longuet-Higgins  except  that 
an  empirically  determined  turbulent  boundary  layer 
solution  was  used  for  the  bottom  matching  condition, 
their  solution  takes  the  form: 


Ifgk.  c»«*  <*.«*♦• 

+  \r  i.a-gE  . 


-0.crr-|-C4^-r]j- 


(2-7) 


where  E«133sinh  kh/a/SD,  /^*</<T/2P,  ^kinematic  viscosity  of 
the  fluid,  and  D  is  the  bottom  roughness. 


For  the  case  of  a  smooth  bottom  the  shape 
function  becomes 


f  ( kh,z/h)*2cosh  2kh  (z/h+1)  + 


3  sin  h  afeA 
*kh 


t (z/h)  -1 ] ( 2-8 ) 


For  large  values  of  kh,  U  — >  0; 
for  small  values  of  kh,  we  have 


U/aVk«0.096/(kh) 


(2-9) 


Again,  for  shallow  water  approximation, 


ittfii 


U//gK=0 .024  *cl 


(2-10) 


Stokes  drift  velocity 

The  classic  expression  of  mass  transport  velocity 
due  to  Stokes  (1847)  is 

U*  <^0sk2kh&+Ll  +  C,  (2-11) 

with  C,  is  an  arbitrary  constant.  For  a  system  with  zero 
total  horizontal  transport  such  as  the  present  case,  we 
must  have 


r  =  _  Sink  iMA 

4-h  s;r?h=kA 


(2-12) 


consequently, 


— iacciA  2kt>($+') 
TshkHd*  I 


TK~  f 


(2-13) 


This  solution  is  also  for  uniform  water  depth.  The 
Stokes  solution  is  for  inviscid  fluid  ana  che  velocity 
profile  is  basically  different  from  that  of  the  previous 
two  cases  as  can  be  seen  in  Fig.  2.10.  the  mean 
offshore  drift  velocity  is 
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u-i/hjSa*  -Jafc‘ 


Ud(z/h) 


(2-14) 


with  the  lower  limit  h(  always  equal  to  zero  and  the 
upper  limit  ht  determined  by  letting  U*0  in  Eq.  (2.13), 
i  .e . , 


h 


*-1  + 


cos A 


(2-15) 


For  the  case  of  shallow  water,  both  terms  on  the 
RHS  of  Eq.  (2.13)  become  large  quantities  and  the 
velocity  distribution  is  subject  to  large  variation  for 
small  numerical  error.  Substituting 

cosh  2kh(z/h+l)=2  sinha kh(z/h+l ) +1 

and  using  a  small  argument  approximation,  Eq.  (2.13)  can 
be  reduced  to 

U»aVk (z/h+1  )a  for  small  kh 

To  satisfy  the  zero  net  drift  condition,  the  profile  must 
be  shifted  such  that 


-fv* 
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U*ai  <r  k  [  ( z/h+l  )4  -1/3} 


(2-16) 


The  mean  offshore  drift  velocity  now  becomes 


U/ai<rk*0.179 


(2-17) 


For  shallow  approximation, 

U/  /gn=0 . 0445Ha  k'1  =1 . 76  ( h/zf  (2-18) 

Numerical  solution  of  stream  functi „n  theory 

Based  on  the  stream  function  theory  by  Dean 
(1965),  a  wave  stream  function  is  defined  as 

aP(x,z)=Lz/T+  2t  X(n)sinh[”~  (h+z)  ]  cos (^~x )  (2-19) 

The  elevation  of  a  particular  streamline,  can  be 

written  as 


z*(x)«TVW/Lf  X(n)  sinh[ ^(h+zt  )  ]  cos (-~x)  (2-20) 
The  particle  velocity  along  the  streamline  is  then 
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u(x,z*)«-jkVX(n)  cosh[^p2(h+z. )  ]  cos(^~x)  (2-21) 

\ 

Averaging  u  for  one  wave  cycle  using  finite  summation, 
the  Lagrangian  mass  transport  velocity  is  obtained 


D(5Jmi  £u,Xi 

i»i  1 


i 


t)At; 


(2-22) 


where  NL  is  the  total  segment  number  of  one  wave  length, 
4t,'»4Xj/Uj,  is  the  time  fraction  the  particle  travelled 
over  Axj  relative  to  wave  motion,  and  ax,-  is  the  length 
of  each  small  segment.  Evaluating  U(: J  for  different 
'Ji,  the  vertical  drift  velocity  profile  can  be 
constructed. 

Details  of  the  ‘Lagrangian  mass  transport  by 
stream  function  theory  are  given  in  Appendix  I.  Applying 
this  numerical  scheme,  the  values  of  hi,  the  upper 
intersect,  and  U  can  also  be  determined.  The  results  for 
the  current  experimental  conditions  are  tabulated  in 
Table  2.2. 


Fig.  2.12  shows  the  dimensionless  depth  h, ,  h* 
for  the  four  theories  discussed  above.  The  data  of  the 
current  experiment  as  well  as  those  of  Russel  and  Osorio, 
and  Mei  et  al.  are  plotted  on  the  same  figure  for 
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Table  2. 

2  Mean 

drift  by 

stream  : 

function 

theory 

case 

(H  /h)b 

V/lgn 

40/ (  Jghk* 

H*)  kh 

U/aVk 

h  z. 

1 

.616 

.0083 

.080 

.645 

.117 

-.380 

2 

.585 

.0120 

.135 

.594 

.213 

-.280 

3 

.771 

.0120 

.222 

.469 

.249 

-.274 

4 

.608 

.0113 

.160 

.531 

.233 

-.233 

5 

.832 

.0115 

.286 

.400 

.239 

-.278 

6 

.814 

.0122 

.243 

.444 

.260 

-.274 

7 

.931 

.0115 

.355 

.360 

.312 

-.278 

8 

.699 

.0129 

.171 

.  550 

.218 

-.270 

9 

.850 

.0127 

.231 

.470 

.246 

-.274 

10 

.910 

.0121 

.340 

.377 

.308 

-.277 

11 

.706 

.0107 

.248 

.415 

.302 

-.288 

12 

.343 

.0098 

.343 

.339 

.336 

-.277 

13 

1.014 

.0118 

.431 

.331 

.341 

-.277 
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Limits  of  offshore  drift  depths,  h^, 
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comparison.  Data  from  sloping  beach  cases  are  found  to 
have  larger  deviation  from  theories  than  those  from 
constant  depth.  Also,  ha  obtained  by  stream  function 
theory  is  not  just  a  function  of  kh  as  suggested  by  the 
other  three  theories.  The  depth  averaged  drift  velocity 
derived  above  are  summarized  in  Table  2.3. 

2.4  Results  and  Comparison 

Since  the  drift  velocity  is  a  second  order 
property,  it  was  normalized  by  aVk.  in  terms  of  shallow 
water  parameters,  the  drift  is  found  to  be  proportional 
to  wave  celerity  and  We  now  proceed  to  compare  these 
results  with  the  experimental  data.  Fig.  2.13  plots  the 
dimensionless  offshcre  drift  velocity,  U/a'c'k,  vs.  kh  on 
a  semi-log  scale.  The  feur  theoretical  solutions  are 
different  in  shape.  Longuet-Higgins'  solution  (LK)  first 
decreases  with  increasing  kh,  reaches  a  minimum  around 
kh«l.l.  For  kh  larger  than  1.1,  the  drift  velocity 
increases  again,  which  is  physically  unlikely  to  happen. 
The  Stokes  solution  (ST)  is  a  monotonically  decreasing 
function  with  increasing  kh  and  the  value  almost 
maintains  constancy  for  kh<l.  The  solution  by  Wang  and 
Liang  (WL)  is  similar  to  Longuet-Higgins  for  kh<l.  For  a 
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larger  value  of  kh,  the  solution  becomes  closer  to  the 
Stokes’  results.  Stream  function  solution  shows  a 
similar  trend  as  that  of  Wang  and  Liang's,  but  generally 
underestimates  U  for  kh<0.5.  However,  it  certainly  is 
better  than  Stokes  theory  at  shallow  water.  Most  of  the 
experimental  data,  as  can  be  seen,  fall  in  between  the 
curves  of  LH  and  WL.  The  WL  curve  seems  to  agree  better 
with  the  data,  at  least  for  kh>0.5.  The  Stokes  solution 
is  apparently  inadequate  for  shallow  water.  In  this 
figure,  the  basic  formulas  ar?  used,  i.e.,  column  2  in 
Table  2.3.  Figure  2.14  compares  the  experimental  results 
with  shalllow  water  approximation,  column  4  of  Table  2.3. 
The  dimensionless  offshore  drift  velocity  is  plotted 
against  It  should  be  noted  here  that  the  parameter 
used  in  the  ordinate  is  different  for  the  Stokes  solution 
as  should  be  evident  from  Table  2.3.  Both  parameters  are 
calculated  for  stream  function  theory  and  are  given  in 
Table  2.2.  The  non-dimensional  parameter  U '/ Jqfi  as 
computed  seems  to  vary  in  a  very  narrow  range  around 
0.012  which  compares  well  with  the  data  range.  Least 
square  fitting  of  the  data  yields  a  slope  of  0.022,  very 
close  to  WL  curve,  which  could  just  be  a  coincidence.  LH 
curve  gives  a  much  higher  estimate  compared  to  the  data. 
Using  the  second  parameter,  \5/(a*  k*  JqE) ,  results  from 
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Table 

Theory 

Longuet-Higgins 
Wang  and  Liang 
Stokes 


2.3  Offshore  mean  drift  velocity 


Basic 

formula 

u 

.  F(feM 

■4  *k  h 

u 

_  pVbh) 

Oi'-Tk 

4  siafkh 

0 

-  F"(feM 

a*<rle 

4-Si^*feh 

Deep 


Shallow 


Z0.0j4kh  -=r  =  0.  060 

m“°24 

U—O 


F(kh)-~  J^{2cosh  2kh  (z/h+l )  +3+kh[ 3  (z/h  f  +4  ( z/h)  +1  ] 

sinh  2kh  +3[sinh  2kh/2kh  +  3/2  ]  [  ( z/h)* -1  ]  }d(  z/h) 
F'(kh)— J^{2cosh  2kh  (z/h  +1)  +32*^  [  ( z/h )a  -1  ] } d ( z/h ) 
F"(kh)*  {2cosh  2kh  (z/h+l)  -  sinh  2kh/kh)d( z/h) 


Fig,2J.3,  Results  of  mean  offshore  drift  velocity  vs.  kh 


Fig.2J.4a  Results  of  mean  offshore  drift  velocity  vs.  k 
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stream  function  show  excellent  agreement  with  the  data. 
Since  the  data  range  is  not  truly  in  shallow  water  range, 
it  is  not  understood  why  this  latter  expression  is  better 
for  stream  function  results. 

Comparing  Figs.  2.13  and  2.14,  one  observes  that 
the  general  formulas  yielded  better  results  than  their 
shallow  water  approximation  except  for  stream  function. 
This  is  partially  explainable  as  the  range  of  kh  tested 
all  fall  in  the  intermediate  water  depth  range  (on  the 
lower  end  close  to  the  shallow  water  case). 

All  the  solutions  presented  here  are  for  constant 
water  depth.  The  effect  of  beach  slop*  is  not  included. 
At  present,  there  is  no  analytical  solution  for  the 
sloping  beach  case  to  facilitate  comparison.  Based  upon 
the  three  data  sets,  offshore  drift  velocity  on  the  mild 
slope  cases,  in  general,  is  higher  than  that  of  the 
steeper  slope  cases.  Fig.  2.15  compares  the  cases  when 
the  drift  velocity  is  plotted  against  kh  to  that  when  the 
drift  velocity  is  plotted  against  a  modified 
dimensionless  parameter,  khs*.  Because  of  the  limited 
cases  (only  three  different  slopes),  conclusions  could 
not  be  made  other  than  that  the  slope  effect  is  probably 
not  as  important  as  the  local  water  depth  and  input  wave 
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conditions . 

2.5  Conclusion 

Laboratory  experiments  were  conducted  to 
determine  the  drift  velocity  at  the  breaking  point  on 
sloping  beaches.  It  was  found  that,  irrespective  of 
breaker  types,  the  drift  velocity  is  in  the  onshore 
direction  near  the  surface  and  close  to  the  bottom.  In 
the  main  flow  column,  the  drift  velocity  is  always 
offshore.  The  offshore  drift  velocity  shows  a  more 
uniform  vertical  distribution  than  that  in  the  offshore 
region  probably  due  to  the  increased  vertical  momentum 
transfer  by  turbulence  near  the  breaking  point.  This 
vertical  distribution  is  very  similar  whether  the  breaker 
is  plunging  or  spilling. 

The  predominant  offshore  drift  velocity  should 
play  an  important  role  in  offshore  sediment  transport. 
Since  the  vertical  distribution  is  rather  uniform,  the 
vertically  averaged  offshore  drift  velocity  is  a 
meaningful  index.  Four  different  constant  depth  wave 
theories  were  used  to  compute  the  mean  vertically 
averaged  offshore  drift  velocity.  The  solutions  of 
Longuet-Higgins  and  Wang  and  Liang  were  found  to  compare 
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reasonably  well  with  experiments  in  the  range  tested, 
with  the  .  latter  being  a  better  fit.  The  Stokes  theory 
performs  poorly.  Results  from  stream  function  are 
significantly  better  than  Stokes'  but  still  underestimate 
the  drift  velocity  in  shallower  water  cases. 

The  test  range  was  rather  limited  both  in  terms 
of  wave  parameter,  kh  and  slope  variations.  Further 
laboratory  studies  should  be  continued  to  cover  wider 
ranges. 

Techniques  for  obtaining  Lagrangian  drift 
velocity  needs  improvement.  The  particle  tracing  method 
used  may  not  produce  enogh  data  at  the  desired  location 
and  generally  produced  very  scattered  results  especially 
near  breaking. 

Reflection  from  beach  was  not  measured,  since  the 
drift  velocity  is  proportional  with  amplitude  squared, 
the  effect  may  be  small  in  comparison  with  the 
uncertainty  in  measurements. 


CHAPTER  3  SUSPENDED  SEDIMENTS 


3.1  Introduction 

Measurement  of  suspended  sediment  in  the  surf 
zone  has  not  been  a  favored  topic  among  investigators. 
Documented  field  data  are  scarce.  7et,  such  information 
is  imperative  to  a  better  understanding  on  the  basic 
mechanics  of  littoral  processes  under  breaking  waves  of 
which  our  present  knowledge  is  rather  rudimentary.  For 
instance,  practically  all  the  existing  longshore  sediment 
transport  formulas  are  based  upon  the  wave  energy  flux 
approach  with  no  real  considerations  of  the  fundamental 
mechanics.  Such  an  approach  leaves  little  room  for 
future  improvement.  As  far  as  onshore/offshore  sediment 
transport  is  concerned,  there  exists  no  predictive 
formula  even  at  the  purely  empirical  level. 

A  number  of  factors  can  be  cited  that  hinder  the 
progress  of  the  research  toward  this  direction.  First  of 
all,  surf  zone  measurement  is  a  tedious  and  difficult 
task.  Secondly,  there  is  no  proven  instrument  to  measure 
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suspended  sediment  in  the  surf  zone;  all  the  existing 
techniques'  appear  to  be  primitive.  Thirdly,  the  most 
useful  data  are  those  collectd  under  extreme  weather  with 
simultaneous  measurement  of  all  the  pertinent 
environmental  factors;  hence,  successful  measurement  is 
often  a  chance  event.  Fourthly,  surf  zone  data  often 
appear  to  be  random  with  large  variance  which  makes 
interpretation  difficult  and  presentation  unimpressive. 

Perhaps  the  earliest  serious  attempt  to  measure 
suspended  sediment  distribution  in  water  column  in  the 
surf  zone  was  by  Matts  (1953).  He  designed  a  pump 
sampler  that  took  time  averaged  samples  by  pumping 
sediment-ladened  water  through  intake  nozzles  situated  at 
various  elevations.  The  sampler  was  laboratory  tested. 
He  stipulated  that  if  the  nozzle  opening  is  of  the  order 
of  1.27  cm  and  the  pumping  velocity  is  approximately 
twice  the  maximum  water  particle  velocity,  the  sample 
collected  should  be  within  about  15  percent  of  the  true 
suspension  on  a  time  averaged  basis.  The  experiments 
were  conducted  from  a  pier  located  at  Pacific  Beach, 
Mission  Bay,  California.  Wave  data  were  partially 
visually  observed  and  partially  measured  by  a  pressure 
gage.  Using  a  similar  device,  Fairchild  (1972)  took  in 
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excess  of  800  suspended  sediment  samples  along  the  City 
Pier,  Ventnor,  New  Jersey,  and  Jennettes  Pier,  Nags  Head, 
North  Carolina.  Most  of  the  data  were  collected  in  the 
surf  zone.  Wave  heights  were  reported  but  the  method  of 
measurement  was  not  mentioned.  Since  both  Watts  and 
Fairchild's  data  were  collected  adjacent  to  piers, 
quentions  have  been  raised  as  to  the  effects  due  to  the 
presence  of  the  structures. 

Pukushima  and  Mizoguchi  (1958),  Hom-ma  and 
Horikawa  (1965)  and  Tanaka  (1975)  tried  suspended 
samplers  made  of  bamboo  poles  to  measure  the  vertical 
distribution  of  suspended  sediment  along  the  Japanese 
coast.  Kana  (1978,  1980),  using  a  water  sampler  device 
similar  to  the  bamboo  sampler,  manually  collected  over 
1500  suspended  sediment  samples  in  the  surf  zone  near 
Price  Inlet,  South  Carolina  and  at  an  experimental  pier 
site  in  North  Carolina.  His  data  should  be  considered 
quasi- instantaneous.  The  wave  information  was  obtained 
by  visual  observations  with  the  aid  of  a  marked  pole. 
Brenninkmeyer  (1974)  used  an  optical  device  termed 
almometer  to  measure  the  sediment  concentration  across  te 
surf  zone  on  the  beach  at  Point  Mugu,  California.  His 
data  were  instantaneous  and  continuous  but  their  absolute 


values  were  not  reliable  because  of  the  uncertainty  in 
calibration.  Morris  (1977)  measured  suspended  sediments 
at  Torrey  Pines  beach,  California  using  a  pump  sampler 
along  with  an  optical  sensor.  The  data  reported  by  these 
investigators  were  either  in  tabulated  form  or  on  log 
papers  to  fit  a  power  or  an  exponential  law  of 
distribution. 

More  recently,  Nielson  and  Green  (1981)  measured 
suspended  sediment  along  the  Australian  coast  using  a 
sediment  sampler  operating  on  hydrostatic  pressure 
difference  instead  of  pumping.  A  more  serious  attempt 
was  m&5e  by  them  to  investigate  into  the  mechanics  of 
sediment  suspension  in  light  of  these  field  data. 

This  chapter  reports  the  results  of  field 
measurements  conducted  along  the  Island  of  Sylt,  Germany. 
Based  upon  the  results  and  those  of  others,  a  new 
approach  is  proposed  here  to  explain  the  mechanics  of 
sediment  suspension  and  to  develop  predictive 
capabilities.  To  augment  the  field  data,  laboratory 
experiments  were  also  conducted  to  further  test  the 
proposed  method  of  computation. 


3.2  Field  Experiment 
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The  field  experiments  described  herein  were  part 
of  a  larger  scale,  multipurpose  exercise.  The  site 
location  and  area  hydrograph  are  shown  in  Pig.  3.1.  The 
site  is  exposed  to  the  North  Sea  and  is  characterised  by 
the  existence  of  an  underwater  transverse  ridge  running 
parallel  to  the  shore  approximately  500  m  from  the 
shoreline.  The  sea  bottom  for  the  general  area  is  mainly 
sand  with  size  distributions  varying  from  fine  to  coarse. 
At  the  test  site,  the  material  is  medium  sand  as  will  be 
shown  later. 

The  beach  slope  varies  from  1:100  under  normal 
wave  conditions  to  as  steep  as  1:50  during  storm  periods. 
The  tide  is  mainly  semidiurnal  with  a  normal  range  of 
about  2  m.  Consequently,  the  beach  front  is  constantly 
assaulted  by  various  wave  conditions  during  one  tidal 
cycle.  In  the  vicinity  and  including  the  test  area,  the 
beach  is  highly  erosional. 

The  first  set  of  data  were  collected  during  and 
right  after  a  storm  on  September  3-7,  1978.  The 
sustained  wind  reached  more  than  25  m/sec  during  the 
height  of  the  storm  and  the  sea  was  fully  risen.  The 
second  set  of  data  were  taken  on  September  26  and  27, 
1979  when  a  moderate  storm  of  lesser  strength  (sustained 
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wind  speed  approximately  11.5  m/sec)  hit  the  area. 

In  both  experiments,  continuous  samples  were 
taken  by  a  pump  sampler  designed  in  accordance  with  the 
specifications  recommended  by  Watts  (1953).  It  has  four 
intake  nozzles,  mounted  on  a  vertical  rack  at  intervals 
of  50  cm.  The  sampler  was  supported  by  a  vertical  staff 
set  in  the  surf  zone.  To  minimize  the  effect  of  the 
supporting  structure,  the  sampler  set  was  set  apart  from 
the  supporting  staff  with  a  horizontal  clearance  of  1  m. 
The  assembly  is  shown  schematically  in  Fig.  3.2.  The 
nozzle  opening  was  2  cm  and  the  pumping  rate  was 
approximately  8.5  liters  per  minute.  samples  were 
collected  at  a  constant  volume  of  36  liters  per  sample  (4 
to  5  minute  sampling  time)  and  were  taken  at  one  hour 
intervals  for  approximately  half  a  tidal  cycle  (6  hours). 
During  the  second  experiment,  the  bottom  nozzle  was  often 
clogged  with  sand  and  these  samples  were  discarded.  In 
the  second  experiment,  in  addition  to  the  contineous 
sampling,  instantaneous  samples  were  taken  by  a  bottle 
sampling  device  modified  from  the  one  first  used  by  Kana 
(1976).  It  consisted  of  a  2  m  long  mounting  pole, 
support  brackets  and  four  2  liter  bottles  closed  off  by 
hinged  doors.  The  whole  assembly  was  mounted  on  a  pole 
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with  sliding  tracks  (Fig.  3.3)  and  the  doors  could  be 
triggered  to  close  remotely  to  enable  instantaneous 
sampling.  Sample  collectors  were  mounted  at  30,  70,  110 
and  150  cm  above  the  bottom.  The  complete  assembly  was 
tested  out  satisfactorily  on  September  26,  1979,  one  day 
ahead  of  the  scheduled  actual  experiment.  Unfortunately, 
an  overnight  storm  broke  the  mounting  pole.  The  original 
plan  had  to  be  abandoned.  Instead,  the  instantaneous 
sampling  had  to  »  carried  out  by  manually  carrying  the 
sampler  into  the  surf  zone  and  triggering  the  device  at 
the  desired  instance.  Because  of  the  high  surf 
condition,  this  manual  operation  was  difficult  to  carry 
out;  the  sampling  location  and  timing  could  not  be 
controlled  as  accurately  as  from  a  stationary  station. 
Despite  this  setback,  a  total  of  sixteen  runs  were  made 
at  time  synchronized  with  continuous  samplings. 

Surface  waves  and  currents  were  measured  at 
locations  shown  in  Fig.  3.4.  The  deep  water  gages  W1 
and  W2  were  bottom  mounted  echo  transreceivers.  These 
echo  sounders  were  capable  of  measuring  water  surface 
variations  to  +/-0.5  cm  for  water  depths  up  to  90  m.  The 
shallow  water  gages  were  staff  mounted  pressure 
transducers  of  type  MDS76  as  manufactured  by  H.  Maihak 
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AG,  Hamburg  with  a  pressure  range  of  0-1  kg/cm  .  The 
current  meters  were  electromagnetic  two-component  type 
manufactured  by  COMEX.  During  the  experiment,  beach  sand 
samples  were  also  collected  and  beach  profiles  were 
surveyed  at  low  tide. 

suspended  sediment  data  collected  during  the 
storm  period  on  Sptember  3,  1978  are  presented  in  Table 
3.1;  the  field  conditions  during  the  same  period  are 
summarized  in  Table  3.2.  Fig.  3.5  displays  the  vertical 
concentration  distribution  for  the  September  3  data. 

The  beach  profile  surveyed  at  low  tide  together 
with  the  beach  sand  grain  size  distributions  at  various 
stations  are  shown  in  Fig.  3.6.  As  can  be  seen  here, 
the  beach  material  was  medium  to  coarse  sand  and  became 
coarser  further  offshore.  At  St.  50,  .he  median  grain 
size  was  about  0.38  mm  which  gradualy  increased  to  0.65 
mm  at  the  offshore  bar.  An  anomaly  was  observed  in  that 
the  material  was  finer  at  St.  120  than  at  St.  110  which 
was  shoreward  of  St.  120. 

The  suspended  sediment  data  collected  by 
continuous  sampler  during  September  26  and  27,  1979  are 
given  in  Table  3.3.  Table  3.4  summarizes  the  results  of 
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Table  3.1  Suspended  sediment  pump  data 
(September  3,  1978) 


NO. 

Cl* 

* 

CM 

U 

C3* 

C4* 

(mg/1) 

(mg/1 ) 

(mg/1) 

(mg/1) 

1 

68 

102 

163 

9 

2 

109 

185 

168 

186 

3 

— 

157 

255 

346 

4 

209 

194 

291 

543 

5 

254 

173 

269 

386 

6 

174 

141 

165 

160 

7 

152 

187 

190 

223 

*  Elevations  of  concentration  measurements  were: 
Cl:  160  cm,  C2:  110  cm,  C3:  60  cm,  C4:  10  cm 
from  sea  bed. 
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Table  3.2 

Field 

condition  (September  3,  1978) 

No. 

Tide 

h* 

H 

(m) 

(m) 

(m) 

1 

0.45 

1.50 

0.66 

2 

0.82 

1.90 

0.98 

3 

1.15 

2.20 

0.95 

4 

1.20 

2.30 

0.93 

5 

1.10 

2.20 

0.93 

6 

1.00 

2.10 

0.87 

7 

0.90 

2.00 

0.69 

*  estimate  at  St.  110  based  on  measurement  at 
low  tide. 


Sep.  3,  1978 
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Table 

3.3  Suspended 

sediment  pump  data 

(Septmber 

26-27, 

1979) 

CASE 

Cl* 

C2* 

C3* 

(mg/1) 

(mg/1) 

(mg/1 ) 

N126 

558 

273 

3940 

N226 

550 

288 

1029 

N327 

747 

971 

841 

N427 

697 

826 

500 

N527 

415 

735 

571 

N627 

553 

644 

332 

*  Elevations  of  measurements: 

Cl:  70  cm,  C2:  120  cm,  C3:  170  cm 
from  sea  bed. 
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Table  3.4  Suspend  sediment  bottle  data 
{September  26-27,  1979) 


CASE 

LOCATION 

Cl* 

C2* 

C3* 

(M) 

(mg/1 ) 

(mg/1 ) 

( mg/1 ) 

1 

120 

941 

225 

2 

120 

479 

270 

280 

3 

105 

240 

74 

96 

4 

120 

530 

200 

5 

105 

390 

263 

6 

105 

117 

7 

90 

189 

126 

93 

3 

90 

221 

102 

9 

90 

293 

70 

10 

90 

64 

11 

90 

766 

134 

149 

12 

90 

641 

86 

33 

13 

90 

127 

348 

14 

90 

240 

15 

90 

449 

73 

65 

16 

90 

561 

161 

* 

Elevations:  Cl 

:  30  cm, 

C2:  70  cm,  C3: 

120  cm 

from  sea  bed. 
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Table  3.5  Field  conditions  (September  26-27,  1979) 


Case 

Tide 

h 

H 

T 

(m) 

(m) 

(m) 

(s) 

N126 

.90 

.85 

.85 

6.0 

N226 

1.06 

.85 

.80 

6.0 

N327 

.91 

.80 

.95 

5.9 

N427 

.99 

.80 

.94 

6.2 

N527 

.95 

.30 

.96 

6.3 

N627 

.39 

.80 

.98 

6.0 

1 

.46 

.65 

2 

.48 

.65 

3 

.77 

.75 

4 

.77 

.65 

5 

.80 

.75 

6 

.81 

.75 

7 

.93 

.30 

8 

.96 

.30 

9 

.57 

.30 

10 

.97 

.30 

11 

.96 

.30 

12 

.94 

.30 

13 

.91 

.30 

14 

.91 

.30 

15 

.89 

.30 

16 

.88 

.30 

j 

J 


I 


J 
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the  instantaneous  suspended  sediment  samples  from 
September  27.  The  field  conditions  during  the  same 
period  are  given  in  Table  3.5.  The  beach  profile  at  low 
tide  together  with  the  beach  sand  grain  size 
distributions  are  shown  in  Fig.  3.7. 

Grain  size  distributions  were  obtained  by  sieve 
analysis  and  by  fall  velocity  analysis.  The  results  from 
both  methods  are  comparable.  The  fall  velocity  analysis 
used  for  this  task  is  connected  to  a  mini-computer  which 
pertains  structured  statistical  analysis  and  plots  size 
distribution.  Detailed  data  tabulation  and  presentation 
can  be  found  in  Wang  et  al,  (1980).  Those  data  revealed 
several  aspects: 

A.  On  the  average,  the  grain  size  of  the  suspended 
sediment  is  much  smaller  than  the  bottom  sand.  For 
example,  the  data  of  September  3,  1978  shows  that  the 
mean  median  grain  size  for  all  the  suspended  samples  is 
0.2  mm  as  compared  to  that  of  0.45  mm  for  the  bottom 
material.  For  the  September  27,  1979  data,  the 
corresponding  values  are  0.2  mm  and  0.6  mm,  respectively. 

B.  The  suspended  sediment  grain  size  distribution  is  much 
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more  uniform  than  the  available  bottom  sand  as  can  be 
seen  from  the  values  of  their  standard  deviation. 

C.  In  the  water  column,  the  sand  sizes  decrease  with  the 
elevation . 

3.3  Laboratory  Experiment 

The  experiment  was  conducted  in  the  same  steel 
wave  flume  deescribed  in  Chapter  2.  Sand  beach  was 
placed  at  the  glass  section.  The  mean  diameter  of  the 
sands  was  0.31  mm  with  fall  velocity  variance  of  0.36. 
The  quantity  of  sand  was  large  enough,  no  overtopping  or 
exposure  of  the  tank  bottom  ever  occurred.  Two 
capacitance  type  wave  probes  were  used  for  wave 
measurements.  One  located  at  the  breaking  point  and  the 
other  at  a  constant  depth  section  5  m  from  the  generator. 
The  constant  depth  probe  was  installed  on  a  moving  cart 
to  facilitate  reflection  measurements.  The  response  time 
for  the  probes  was  about  0.03  seconds. 

Horizontal  and  vertical  particle  velocities  at 
the  breaking  were  measured  using  a  Marsh-McBirney 
electromagnetic  current  meter,  model  523M  with  0.1  second 
time  constant.  The  sensing  volume  was  a  sphere  of  1.5" 


diameter  (three  times  of  the  sensor  diameter);  data  so 
collected  were  thus  an  average  over  this  distance.  The 
sediment  concentration  was  measured  by  _,n  Iowa  Sediment 
Concentration  Measurement  System  (ISCMS)  designed  at  the 
University  of  Iowa,  the  major  problem  of  this  instrument 
is  the  intermittent  signal  output  as  widely  discussed  in 
the  literature  (Nakato  et  al.,  1977;  Glover,  1977; 
Locher  et  al.,  1976).  The  sensing  volume  is  about  1  mm 
sphere  and  the  frequency  response  is  1  KHz. 

The  data  were  fed  into  an  on-line  Digital  MINCH 
minicomputer  through  A/D  converter.  The  sampling  rate 
and  the  sample  size  were  fixed  at  64  Hz  and  64  seconds, 
respectively.  Wave  record  at  constant  depth  was 
registered  by  a  HP  7402A  chart  recorder.  A  Tektronix 
Oscilloscope  with  model  3A74  4-trace  amplifier  was  used 
for  monitoring  the  A/D  data  collection. 

At  the  beginning  of  the  experiment,  a  constant 
slope  sand  beach  was  constructed  and  this  beach  was  then 
subject  to  wave  action  for  1C  to  12  hours  until  a 
quasi-equilibrium  condition  was  reached.  At  the  breaking 
point,  simultaneous  data  of  surface  elevation,  horizontal 
and  vertical  velocities  and  concentrations  were  then 
taken.  The  measurements  of  the  last  three  quantities 


69 


started  from  the  lowest  elevation  which  was  just  above 
the  ripple  crest,  then  progressing  upward  at  0.5  cm 
interval  for  the  first  3  levels,  and  1  to  2  cm  intervals 
afterward  depending  on  the  local  water  depth.  A  total  of 
5  to  9  elevations  for  each  case  were  taken.  Table  3.6 
summarizes  the  test  conditions  of  the  experiments. 

3.4  On  the  Theory  of  Vertical  Concentration  Distribution 

Of  great  interest  to  sediment  transport  is  the 
vertical  distribution  of  suspended  sediment  in  the  water 
column.  The  problem  of  sediment  suspension  in 
oscillatory  flow  is  an  unsettled  issue.  A  general 
examination  on  the  mechanics  of  sediment  suspension  in  a 
wave  field  is  given  first  before  the  analysis  cf  the 
actual  field  and  experimental  data. 

At  present,  sediment  suspension  in  a  fluid  media 
is  treated  as  a  diffusion-dispersion  process,  and  the 
governing  equation  takes  the  following  form: 

SC/dt+V.  (Cv^)-7.  (PVC)  (3-1) 

where  C  is  the  volumetric  or  weighted  sediment 
concentration;  v^  is  the  particle  velocity  vector,  *•'  is 


Table  3.6  Wave  conditions  of  sediment  experiment 


CASE 

T  He  hc 

Hk  h.  kh 

R+  I04§ 

Ur 

\r 

(s)  (cm)  (cm)  (cm)  (cm) 

<%> 

(cm) 

(cm) 

A-04 

2.05  6.035. 

0  12.5  12.5  .37 

6 

1 2 1 

1.53 

2.0 

15.0 

A-05 

2.69  8.3  34. 

5  13.7  17.0  .31 

21 

93 

1.91 

4.0 

15.0 

A-06 

2.10  7. 7  36. 

8  12.0  15.8  .38 

10 

114 

1.28 

1.0 

10.0 

B-06 

1.79  5.2  36. 

8  10.8  13.8  .43 

26 

136 

1.08 

2.0 

10.0 

A-07 

1.59  6.3  36. 

8  11.0  13.8  .49 

10 

161 

0.88 

2.0 

8.0 

3-07 

1.43  6.3  36. 

8  10.0  15.0  .57 

12 

157 

0.52 

2.0 

8.5 

A-08 

1.23  8.4  36. 

8  10.0  17.0  .73 

5 

203 

0.33 

1.5 

7.0 

B-OS 

1.08  11.1  36. 

8  11.0  21.5  .99 

■> 

247 

0.16 

1.0 

6.0 

C-08 

1.08  11.1  36. 

8  10.2  9.0  .59 

•> 

344 

1.06 

7.0 

50.0 

A-09 

2.21  5.1  35. 

7  10.0  16.5  .38 

9 

88 

1.02 

1.0 

10.0 

A-10 

3.74  9.0  34. 

0  17.0  19.0  .24 

20 

66 

3.13 

1.0 

10.0 

B-10 

3.07  10.9  34. 

5  16.5  18.5  .29 

6 

76 

2.11 

1.0 

10.0 

A-ll 

3.44  6.5  34. 

0  14.0  18.0  .25 

? 

66 

2.64 

2.0 

15.0 

(Table  3.6  continued, 

el  s' 

nations ) 

case 

zl  z2 

z  3  z4  z5 

z6 

zl 

zB 

z9 

A-04 

-.97  -.94  - 

.87  -.86  -.81  - 

.41 

A-05 

-.96  -.93  - 

.94  -.92  -.80  - 

.62 

-.35 

A-06 

-.95  -.92  - 

.85  -.79  -.69  - 

.29 

B-06 

-.96  -.93  - 

.89  -.78  -.67  - 

.33 

A-07 

-.96  -.93  - 

.91  -.82  -.71  - 

.29 

B-07 

-.97  -.93  - 

.93  -.83  -.73  - 

.60 

-.47 

-.33 

A-03 

-.91  -.94  - 

.38  -.82  -.74  - 

.62 

-.50 

-.25 

-.37 

B-08 

-.95'  -.93  - 

.38  -.79  -.65  - 

.47 

-.23 

C-08 

-.94  -.89  - 

.83  -.72  -.67 

A-09 

-.97  -.94  - 

.88  -.76  -.64  - 

.49 

A-10 

-.97  -.95  - 

.39  -.79  -.69  - 

.53 

-.38 

-.27 

3-10 

-.95  -.95  - 

.89  -.84  -.82  - 

.65 

-.49 

-.32 

A-ll 

-.99  -.95  - 

.92  -.81  -.64  - 

.42 

+  reflection  coefficient  from  constant  depth  measurement 
?  not  measured 
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the  molecular  diffusion  coefficient.  For  an  oscillatory 
flow  field,  the  concentration  and  velocity  can  both  be 
divided  into  three  components: 


C=C+CW+C' 


(3-2) 


and 


v.j=v+vw+v’  (3-3) 

Here  C  and  v  are  the  mean  values,  Cw  and  vw  are  the 
oscillatory  components,  and  C',  v'  are  the  turbulent 
fluctuations.  Substituting  Egs.  (3.2)  and  (3.3)  into 
Eq.  (3.1),  and  taking  time  average  yields 

7.(Cv+Cwvv+C'v')^.(WC)  (3-4) 

If  the  molecular  diffusion  is  small  and  the 
gradients  in  the  horizontal  plane  are  much  smaller  than 
those  in  the  vertical  direction,  Eq.  (3-4)  becomes 


(3-5) 
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here  z  is  the  vertical  axis,  with  positive  pointing 
upward,  v,  v„  and  v'  are  the  vertical  components  ofj^r,  vC„ 
and  v',  respectively.  Since  there  is  no  net  flux  itiross 

P 

the  free  surface,  Eq.  (3-5)  when  integrated  once  becomes 


Cv+Cwv^+C ' v '  =0 


(3-6) 


Here  the  mean  velocity,  v,  is  usually  negative  and  is  the 
mean  settling  velocity  of  sediment  in  calm  water.  In  the 
usual  context  of  turbulent  diffusion  processes,  the  last 
term  in  Eq.  (3-6)  is  assumed  to  be 


C'v'*- £ 


(3-7) 


where  £  is  the  turbulent  diffusion  coefficient, 
(3-7)  now  becomes 


Cv+CL  vw  -  £  -~- 


(3-8) 


I f  we  let  w6  be  the  magnitude  of  particle  fall  velocity 
in  calm  water,  Eq.  (3-8)  becomes 


-Cw#  +Cwv  -£  0 


(3-9) 
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The  first  and  last  terms  in  Eq.  (3-9)  constitute  the 
conventional  one-dimensional  turbulent  diffusion 
equation;  the  middle  term  represents  the  concentration 
being  convected  by  the  oscillatory  velocity  component  in 
the  vertical  direction.  If  the  terms  of  C„ vw  and  £  are 
known  or  given,  Eq.  (3-9)  can  then  be  used  to  solve  the 
vertical  sediment  distribution.  Unfortunately,  none  of 
these  terms  are  known  in  oscillatory  flow  field;  various 
assumptions  have  to  be  made  to  arrive  at  distribution 
laws  that  are  to  be  verified  by  experiments.  A  number  of 
cases  that  have  been  proposed  by  various  investigators 
are  examined  here. 

Up-to-date  research  efforts  have  been  mainly  on 
determining  £  in  wave  field  neglecting  the  contribution 
due  to  C«Vw  (at  least  not  explicitly  treated  as  a  special 
term).  Then,  depending  upon  the  assumptions  made  on  £, 
various  distribution  functions  can  be  obtained.  Most 
investigators,  for  instance,  treated  both  w0  and  i  as 
constants  which  results  in  an  exponential  distribution. 
Laboratory  experiments,  notably  Bhattacharya  (1971), 
seemed  to  suggest  a  power  law  distribution  which  implies 
that  t  varies  linearly  with  the  vertical  coordinate  z. 
Hom-ma  and  Horikawa  (1963),  following  von  Harman's  mixing 
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length  hypothesis,  derived  an  expression  for  to  vary  as 
sinhJ  kz/sinhJ  kh  which  is  approxitnatly  zJ  in  shallow 
water;  the  resulting  distribution  then  becomes  a 
log-exponential  curve  of  the  form  In  C*C, z *.  Other 
suggestions  were  also  made  by  different  investigators 
(Dally,  1980;  Yang,  1981;  Nielson,  1981,  for  instance) 
to  arrive  at  different  distributions.  Table  3.7  gives  a 
summary  of  various  formulas  on  sediment  vertical 
distributions  as  being  proposed. 

All  the  formulas  given  in  Table  3.7  will  yield 
distribution  functions  that  are  monotonically  decreasing 
with  elevation,  i.e.,  sediment  concentration  decreases 
monotonically  from  the  bottom  to  the  surface.  A  number 
of  experimental  results,  including  some  cases  reported 
here,  revealed  such  not  to  be  the  case.  The 
concentration  can  actually  be  higher  at  surface  layer 
than  mid-depth.  There  are  three  possible  causes  for  such 
anomaly: 

1.  Oscillatory  velocity  field  resulting  in  a  reduced 
effective  fall  velocity  for.  sediment  particles  (Fig. 
3.8a), 


Sources 


Table  3.7  Summary  of  £  and  C  prof i lees 
£  C*(z)=C(z)/Co 

constant  exp[w(z-Za)/£] 


Kz  (  +  ) 
Khu* 
Kvwl 


( z/z  0)** { ~w/K) 

exp[ (-Kwexp[ (z-z0)/u*h] 


ttanh  |£/tanh  ^  3 


k* 


exp[-  f  (f  iJntfhfe  ] 


Dally,  1980 

Wang  and 

Liang,  1975 

Hom-ma  and 


Horikawa,  1962 


£b+£3(z/h)1 


exp[ 


— w  K 


tar T'i 


lS 


_2_ 

k 


Nielsen  and 
Green,  1981 


Ku 


exp[ 


-  W  JK'hKO 

f 3a  <T 


ft**  £ 


-  I  ~ 


*  K:  a  constant 

u  ,  v  :  horizontal  and  vertical  particle  velocities, 
respectively 

u*:  frictional  velocity  («>[T/p ) 
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Pig.  3.8  Conceptual  sketch  of  mechanisms  that  gernerate 
high  concentration  at  higher  elevation 
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2.  Large  horizontal  sediment  concentration  gradient 
coupled  with  favorable  mean  flow  profile  such  as 
illustrated  in  Fig.  3.8b, 

3.  The  existence  of  stationary  eddies  or  circulations  in 
the  surf  zone  to  transport  high  sediment  concentration 
water  to  the  surface  such  as  illustrateed  in  Fig.  3.8c. 

Certainly  the  combination  of  these  causes  could 
also  yield  higher  sediment  concentration  at  the  surface 
layer.  In  this  section,  the  effect  of  Q5„on  the  vertical 
concentration  distribution  is  pursued  further. 

Settling  velocity  of  sediment  particles  in 
oscillatory  flow  field  has  been  found  to  be  smaller  than 
that  in  calm  water  by  a  number  of  investigators.  Ho 
(1964)  showed  that  under  large  acceleration,  fall 
velocity  can  be  reduced  by  as  much  as  70  percent. 
Murray's  result  (1970)  showed  that  under  a  normal 
turbulent  river  condition,  40  percent  reduction  in  fall 
velocity  is  possible.  Wang,  S.-Y.  (1981)  and  Lee  (1981) 
conducted  experiments  on  particle  fall  velocity  in  wave 
field.  Their  results  yielded  large  reduction  in  fall 
velocity  near  the  surface;  this  reduction  diminishes 
exponentially  as  particles  settled  further  downward. 


a 


Wang,  S.-Y.  (1981)  introduced  a  so-called  effective  fall 
velocity  concept  and  was  able  to  produce  vertical 
concentration  profile  that  has  a  minimum  at  mid-depth. 
Nielsen  (1979)  also  discussed  the  fall  velocity  reduction 
in  a  fluctuating  flow  field  but  did  not  pursue  further  in 
his  analysis.  If  Cwv,,  term  represents  a  perturbation  on 
Cw  and  the  presence  of  nonuniform  vertical  velocity 
field  is  to  cause  reduction  in  particle  fall  velocity,  we 
assume  the  following  relationship: 


-Cw  +C  v  =-Cw  (3-10) 

O  "V  Vv 

where  w  is  the  effective  fall  velocity  defined  as 

W3W0 -R  W  (3-11) 

where  R  is  a  fall  velocity  reduction  coefficient  and  w  is 
the  amplitude  of  the  nonuniform  vertical  flow  velocity. 
In  a  wave  field,  we  have  W»v„a  "amplitude  of  vw. 

Substituting  Eq.  (3-10)  into  Eq.  (3-9),  we  have 


He 


-Cw 


[ii 
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This  equation  has  the  same  form  as  the  conventional 
vertical  diffusion  equation  except  that  w  is  no  longer  a 
constant  but  varies  with  vw.  Integrating  this  equation 
yields  the  following  concentration  relationship: 

C(z)  =  C0  £xp[-J*  £(w0-RW)dz]  (3-13) 

where  Ce  is  the  sediment  concentration  at  reference  level 
zQ .  For  convenience,  the  overbar  is  dropped  from  this 
point  on.  The  integration  can  be  carried  out  if  w« ,  R 
and  L  are  given.  Furthermore  if  both  CD  and  ze  can  be 
determined,  C(z)  is  then  quantatively  solved.  The 
problems  of  selecting  £,  R  and  z©  are  dealt  with  in  the 
following  sections. 

3.5  Diffusion  Coefficient  Selection 

As  indicated  in  Table  3.7,  there  are  a  great 
number  of  different  assumptions  on  £..  After  examining 
all  possibilities,  the  author  concluded  that  a  modified 
hypothesis  of  Wang  and  Liang  (1975)  appears  to  be  most 
plausible.  Wang  and  Liang's  original  assumption  is  that 
the  diffusion  coefficient,  £,  is  proportional  to  a  length 
scale,  1,  which  they  treated  as  constant,  and  a  velocity 
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scale,  v,  which  they  assumed  to  be  proportional  to  the 
vertical  component  of  the  wave  motion.  Here,  the 
modified  assumption  uses  the  same  velocity  scale  but  the 
length  scale  is  no  longer  treated  to  be  constant  but 
proportional  to  the  vertical  excursion  of  the  water 
particle  motion,  i.e., 

£  *klv  *k  vWo(z)  *j(z)  (3-14) 

where  \(z)  is  the  amplitude  of  the  vertical  excursion  of 
the  water  particle  and  k  *s  a  constant  of 
proportionality.  For  linear  wave  field,  we  have 

£  =k.HJT sinh2  kz/( 2  sinh1kh)  (2-15) 

In  the  breaking  zone  the  eddy  size  is  likely  to  be 
governed  by  the  local  wave  height  or  local  water  depth, 
therefore,  the  suitable  representation  for  at  breaking 
zone  or  area  of  strong  agitation  is 

6  »K.H*-crsinh  kz/(2  sinh  kh)  (3-16a) 

or 


I  ■fcHhO'sinh  kz/(2  sinh  kh) 


(3-16b) 


In  the  following  paragraphes,  the  concentration  profiles 
under  strong  agitation  {breaking  and  surf  zone)  and 
normal  agitation  will  be  developed  in  parallel  using  Eqs . 
(3-15)  and  (3-16a). 

Eq.  (3-13)  can  now  be  integrated  to  yield  the 
concentration  profiles: 


C/Co* ( 


-hu,*,  ' 


Rsi*4ikk 

K-Wfc 


Exp[ 


2.  Wo.S'wIV'fek 


hj  cr  fe 


(cth  kz-cth 


kz0  )  ] 
(3-17a) 


for  nonbreaking ,  and 

[j  g  2  Ufc  fanA  kh 

C/Co*Exp[R  (2-Z0)/<iH](~^-g-)  K'HlTte  (  3-1 7b ) 

for  breaking  (and  surf  zone). 

This  solution  is  similar  to  Wang  and  Liang's 
except  that  the  first  part  on  the  right  hand  side  which 
is  due  to  the  modification  of  the  effective  fall 
velocity.  For  shallow  water,  Eq.  (3-17b)  reduces  to: 

2_ujo_£> 

C/Co*Exp[R  (z-z,  )/k<H](-§£) 


(3-18) 
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In  the  above  integration,  we  is  assumed  to  be  a  constant, 
which  holds  only  for  uniform  sediment  sice.  In  nature, 
beach  sands  are  nonuniform.  Consequently,  w4  cannot  be 
treated  as  a  constant.  Nielsen  (1979)  considered  this 
problem  by  assuming  that  the  fraction  of  suspension  of 
each  sand  component,  with  corresponding  fall  velocity 
wt  ,  is  proportional  to  the  percentage  of  the  component  in 
the  bed  material,  and  each  component  follows  the  same 
dispersion  law,  i.e., 

C^(2,wt)  =  (w-  )Exp[-J^o  -g-  (w<  -rtW)dz  j  (3-19) 

where  C^.( wx-  )  is  the  fraction  of  the  ith  sand  component  at 
reference  level  z0.  Introducing  the  following 
dimensionless  par?  «  .ers, 

w'*wv/w;  (3-20) 

and 


C (W  )  «  C.^  (V  )/C0  (3-21) 

where  Cc  is  the  total  sediment  concentration  at  the 
reference  level  z,  and  w"0  is  the  mean  fall  velocity.  Let 
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£(w’)dw'  represnts  the  probability  density  function  of 
Cai/  Eq.  (3-19)  can  then  be  integrated  to  obtain  the 
concentration  distribution  as  follows: 

C/Cc  ■£*  f  (W  )Exp[-/^-^(w0  w'-RW)dz]dw’  (3-22) 

If  €  and  R  are  assumed  to  be  independent  of  w’ ,  we  have 


C/Co«  f  (w’ )Exp[-^w’-(33dw’ 


(3-23) 


where  oi  w^  dz  and  (3-J^-^RWdz.  A  number  of 
investigators  have  suggested  a  T-distr ibution  function 
for  the  sediment  distribution,  such  that 


(  I/V  )'X'  w 

f  (w’  )-— -—^-ExpI-w'/V]  v  for  0<w ’  <  (3-24) 

rc  w  ) 


with  V-the  variance  of  the  distribution  function. 

Substituting  Ecs.  (3-24)  and  (3-16a)  or  (3-15) 
into  Eq.  (3-23),  the  explicit  solutions  are  obtained  for 
sediment  with  nonuniform  size  distribution: 


C/C0.(^^-) 


[l-v 


2  W.  S^kfy  (<ed  ki-Ct+A  krb  ) 


te  H* <?*■ 


(3-25a) 


"•!  9a .... 
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for  nonbreaking ,  and 


C/Co=Exp[H  (z-za  )/KiH]  [  1  +  - 


'5.  I ink  kh  '/v 


k,  H'J  ftH  4-fiviti 


10  4 -r'i  ^ 


(3-25b) 


for  breaking.  In  shallow  water.  Eg.  (3-25b)  reduces  to 


C/Co=Ex?[ R  (^J/tK.H/hJJtl^-^-^ln^] 


2  y  _  /  w, , 


rkV 


(3-26) 


In  Eq.  (3-26),  we  identify  the  following  dimensionless 
exprssion: 

C/Co=f  ( z/h ,  z o/h ,  H/h,  w^VHtr,  V)  (3-27) 

Under  breaking  condition,  H/h  fluctuates  in  a  narrow 
range,  the  dimensionless  groups  can  be  further  reduced  to 

C/Co* f( z/h,  z<>/h,  w0  /H O’,  V)  (3-28) 

Fig.  3.9  illustrate  respectively,  the  effects  of  Zo/h, 
vTo/HO-  and  V  on  the  vertical  sediment  distribution.  It  is 
clear  that  z0,  as  a  mathematical  reference  level,  does 


not  affect  the  concentration  profile  shapes,  parameters 
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wo/H(r  denotes  the  agitation  level;  the  higher  agitation 
(lower  w"e/H<r),  the  more  uniform  distribution.  V,  the 
sediment  distribution  variance,  denotes  the  availability 
of  the  suspendable  sediments;  the  poorer  the  sorting 
(larger  variance),  the  more  abundant  the  light  material 
that  could  be  carried  to  a  higher  elevation,  therefore 
showed  more  uniform  distribution  beyond  bedload  level. 

3.6  Fall  Velocity  Reduction 

One  of  the  important  assumptions  made  in  the 
derivation  of  the  vertical  sediment  concentration  profile 
is  that  the  mean  sediment  particle  fa’l  velocity  is 
reduced  due  to  wave  agitation.  The  magnitude  of 
reduction  is  assumed  to  be  proportional  to  the  amplitude 
of  the  vertical  water  particle  velocity,  i.e., 

w/w0  *  1  -  RW/we  (  3 ~ 2 9 ) 

The  value  of  R  must  first  be  established.  A  simple 
analysis  based  upon  consideration  of  change  in  drag  force 
in  oscillatory  flow  is  presented  here. 


I 

I 
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Under  equilibrium  condition,  the  force  balance  on 
a  particle  in  calm  water  can  be  expressed  as: 

G  =  CpA  w//2  (3-30) 

where  G  is  the  gravitational  force;  C-p  is  the  drag 
coefficient;  A  is  the  projected  surface  area  and  we  is 
the  fall  velocity  in  calm  water.  In  oscillatory  flow, 
the  above  equation  is  modified  as 

a^jfq;(w-W  sin  t)lw-W  sin  tl  dt  *  G  (3-31) 

where  is  the  instantaneous  drag  coefficient;  w  is  the 
fail  velocity  in  oscillatory  flow  and  W  is  the  amplitude 
of  vertical  velocity  component.  Equating  Eqs.  (3-30) 
and  (3-31)  and  solving  the  w/w0 ,  we  obtain 

w/w0-[~r  -  0.5  (W/w.  )*  (3-32) 

M> 

for  W/w  <  1,  and 

w/w«-  -~*[l+0.5(W/w)*sin~'(W/w)+1.5  /(w/w-i)f  '*(3-33) 
for  W/w  >  1. 
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Thus,  the  reduction  in  fall  velocity  in 
oscillatory  flow  is  a  function  of  Co/ and  W/wa  .  Fig. 
3.10  plots  this  functional  relationship  with  various 
Cp/CT  ratios  together  with  the  available  laboratory  data. 
The  data  showed  a  large  scatter.  It  is  difficult  to 
claim  any  of  these  curves  is  the  correct  one.  On  the 
other  hand,  it  is  evident  that  the  fall  velocity 
reduction  in  oscillatory  flow  is  quite  significant, 
particularly  when  W/wa  becomes  large.  Consequently,  in  a 
high  agitation  field,  the  reduction  in  fall  velocity  is 
expectd  to  be  larger  than  most  of  the  data  presented. 
For  the  data  reported  in  the  study,  the  values  of  W/w, 
were  in  the  order  of  15  for  the  laboratory  conditions  and 
even  larger  (of  the  order  of  30-40)  for  field  conditions. 
To  simplify  our  computations,  we  treated  R  as  a  constant 
and  set  it  to  be  equal  to  0.07  for  W/w, <14. 3  by  data 
fitting.  For  extremely  high  agitaiton,  the  sediment 
particle  becomes  passive  and  will  be  treated  as  neutral 
buoyant  beyond  W/w0>14.3.  This  resulted  in  a  straight 
line  relationship  between  w/we  and  W/w0  such  as  shown  in 
Fig.  3.10.  If  this  cut-off  velocity  of  W/w0-14.3  is  not 
made,  negative  fall  velocity  occurs,  which  results  in  the 
inversion  of  the  concentration  profiles  beyond  this 
critical  velocity  (see  Fig.  3.15,  Test  12  and  13).  The 


concept  of  negative  fall  velocity  needs  more  careful 

examination,  however. 

3.7  Determination  of  Reference  Level  z0 

The  vertical  sediment  concentration  profile 
discussed  above  is  based  upon  diffusion  consideration  and 
is,  therefore,  only  valid  for  suspended  load  in  the 

region  above  the  bed  load  layer.  The  magnitude  of  the 
suspended  sediment  concentration  (Fig.  3.9)  depends  upon 
the  value  z0  and  the  corresponding  concentration  CG . 
Mathematically,  zQ  can  be  chosen  arbitrarily  in  the 

suspended  sediment  region  as  long  as  Z0  can  be 

determined.  Physically,  this  implies  z„  should  be  chosen 
at  the  boundary  where  C,  can  be  given  as  a  boundary  value 
or  determined  by  considerations  other  than  suspended  load 
equation.  In  river  flow,  Einstein  (1950)  suggested  z„  be 
chosen  at  the  water  edge  of  bed  load  layer.  For  sand 
bed,  he  suggested  z0«2d,  with  d*  diameter  of  the  sand 
particle.  This  criterion  runs  into  difficulty  for  fine 
sand  or  silt  material;  it  also  becomes  impractical  for 
rough  or  rippled  bed-  Nielsen  (1979),  Nielsen  and  Green 
(1981)  among  others  chose  the  ripple  crest  level  for 
reference  which  may  be  inside  the  bedload  layer. 
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Besides,  measurement  of  concentration  at  that  level  was 
generally  scarce  and  had  to  be  extrapolated,  subject  to 
personal  interpretation  and  large  error  due  to  the  scarce 
measurements  at  the  vicinity.  In  this  study,  the  author 
chose  the  reference  level  at  the  boundary  layer  elevation 
using  equations  of  Jonssen  and  Carlsen  (1976): 

S/§b-0.072(kr/|b),/4  (3-34) 

here  f b  is  the  amplitude  of  the  horizontal  particle 
excursion  at  bottom  and  kr  is  the  Nikuradse  roughness 
taken  as  an  arbitrary  number  as  120  d  followed  Nielsen 
and  Green.  For  rippled  bed,  the  following  modified 
equation  for  kr  was  suggested  by  them: 

kr«  120  d  +  1.27  %  (3-35) 
with  1r*ripple  height. 

Based  upon  this  reference  level,  the  reference 
concentration  data  obtained  by  various  investigators  are 
plotted  against  Shields  parameter  in  Fig.  3.11.  The 
Shields  parameter  is  defined  as 


Cot  9^1  ) 


'V  «T/jOgd(s-l) 


(3-36) 


where 

7  =  i-fwU^/2  (3-37) 

a  (44 

fw*Exp[5.213(kr/£b  )  •  -5.997]  for  Vkr>1.7 

*0.28  for  £/kr<1.7  (3-38) 

by  Jonsson  (1966),  uk  is  the  horizontal  particle  velocity 
of  wave  motion  at  bottom.  The  data  shown  here  are  rather 
scattered;  thus  quantitative  dependency  on  ^  is  not 
clear.  Note  that  the  scatter  of  C,  - f  relationshop  did 
not  improve  when  only  the  surface  drag  were  considered 
for  the  Shields  parameter  calculation,  for  example,  see 
Nielsen,  1979  and  Nielsen  and  Green,  1981. 

3.8  Field  and  Laboratory  Data  Presentation 

The  analysis  presented  in  the  previous  section 
presumes  that  when  the  fluid  is  in  a  state  of  high 
agitation  such  that  the  flow  is  fully  turbulent,  the 
characteristic  eddy  scale  is  proportional  to  the  wave 
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height;  while  in  a  less  agitating  state,  the  eddy  size 
is  proportional  to  the  amplitude  of  vertical  particle 
motion.  This  study  emphasizes  the  former  condition, 
field  data  collected  in  the  surf  zone  and  laboratory  data 
collected  at  wave  breaking  point  and  surf  zone  were  used. 
The  conditions  under  which  these  data  were  obtained  have 
been  described  in  the  previous  sections.  In  addition,  a 
few  sets  of  data  from  Nielsen  and  Green  were  also 
included. 

Thirteen  sets  of  laboratory  data  are  given  in 
Fig.  3.12,  together  with  the  theoretical  curves  in 
accordance  with  Eg.  (3.17b).  The  value  is  data 
fitted;  two  values  of  0.4  and  0.25  are  given  to  show  the 
range  of  variations.  The  correlation  coeffient  of  the 
theoretical  values  and  the  measured  quantities  is  used  as 
a  measure  of  the  goodness  of  fit.  The  correlation 
coeffient  is  defined  as 

Y-Cov ( C-r , CM) /YVar ( CT )  .Var(CH)  (3-39) 


where  Cov(CT,CM)  is  the  covariance  of  functions  CT  and 
C^f  and  Var(q)  is  the  variance  of  function  q.  CT  is  the 
theoretical  results  and  CM  is  from  measurement. 
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Fig.  3.13  shows  the  distribution  of  the  goodness 
of  fit.  It  is  seen  that  10  out  of  13  cases  show  .80  or 
higher  correlation.  Fig.  3.14  presents  the  results  of 
the  North  Sea  data.  Since  no  measurements  at  close  to 
the  bottom  was  made,  reference  concentration  evaluation 
is  difficult.  An  uniform  1.5  gram/liter  is  assigned 
arbitrarily  to  facilitate  computation.  Although 
quantitative  conclusion  is  difficult  to  draw  with  only 
three  data  points  in  each  profile,  the  trend  does  show 
the  high  agitation,  large  sand  variance  type  distribution 
(Fig.  3.9b,  c)  as  characterized  at  the  measuring  site. 

A  total  of  21  sets  of  data  from  Nielsen  and  Green 
were  computed.  The  correlation  distribution  is  also 
shown  in  Fig.  3.13.  17  sets  of  them  are  above  0.8. 
Fig.  3.15  shows  some  of  the  typical  results. 

In  the  last  section,  Eq.  (3-17a)  is  suggested 
for  the  concentration  computation  for  weak  agitation 
conditions.  Although  thorough  testing  was  not  done,  a 
brief  computation  of  the  ratio  between  nonbreaking  and 
breaking  concentration  under  similar  wave  condition  is 
compared  to  the  measurements  of  Nielsen  1979,  Fig.  6.25; 
reproduced  as  Fig.  3.16  here).  Since  wave  condition  was 
not  complete  in  the  original  report,  the  following 
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Fig.  3.13  Goodness  of  fit  of  the 
experimental  data  with 
the  theoretical  results 
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approximation  is  used:  kh=0.3,  vTo/H<r*0.1,  H/h*0.5, 
kH*0.02,  zo/h*0.025,  V=0.3.  The  result  shows  a  ratio  of 
1/100  at  mid-depth  which  is  very  close  to  what  was 
measured. 

Fig.  3.17  shows  the  comparison  of  the  present 
model  with  the  models  developed  by  Nielsen  and  Green, 
Wang  and  Liang,  Hom-ma  and  Horikawa,  and  Dally.  It  can 
be  seen  that  for  the  high  agitation  cases,  only  the 
Nielsen  and  Green  model  is  comparable  with  the  present 
results.  Their  model,  however,  introduced  a  lot  of  data 
fitted  cofficints  with  large  uncertainties  while  the 
present  model  used  least  number  of  data  fitted 
coefficients. 

As  a  concluding  remark  to  this  chapter,  sediment 
suspension  under  wave  agitation  can  be  properly  modelled 
as  diffusion  processes.  The  turbulent  viscosity  is 
closely  related  to  the  wave  properties  and  degree  of 
agitation.  The  suitable  velocity  scale  is  the  vertical 
water  particle  velocity,  the  length  scale  is  the  wave 
height,  or  water  depth  for  shallow  case;  under  strong 
agitation,  the  length  scale  is  uniform  over  the  whole 
depth,  while  at  weaker  disturbances,  the  scale  decays 
with  water  depth,  therefore,  the  amplitude  of  the 
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vertical  excursion  is  more  appropriate.  Sediment  grading 
as  well  as  strength  of  agitation  affects  the 
concentration  profiles,  the  stronger  agitation  and  the 
poorer  the  sorting,  the  more  uniform  distribution  over 
the  depth  occurs.  The  reference  level,  however,  only 
affects  the  absolute  values  and  any  level  that  will  give 
accurate  measurement  is  appropriate.  The  effective  fall 
velocity  under  oscillatory  flow  modifies  the  distribution 
and  needs  to  be  properly  considered.  The  present  linear 
model  is  only  a  first  order  approximation,  further 
studies  are  desired  for  a  better  understanding  of  this 
problem.  The  relationship  of  the  reference  concentration 
to  the  flow  and  sediment  properties  is  still  uncertain 
and  need  more  detailed  investigation. 


CHAPTER  4  WAVE  AND  VELOCITY  SPECTRA 

4.1  Introduction 

In  this  study,  all  the  waves  were 
monochromatically  generated.  At  breaking,  however,  the 
spectra  of  these  waves  are  broadened  by  nonlinear 
interaction.  The  spectra  do  not  show  a  fixed  slopes 
which  is  typical  for  waves  in  finite  water  depth.  Huang 
et  al.  (1981)  developed  a  spectral  model  for  general  sea 
state  in  deep  water.  The  spectral  slope  is  a  function  of 
Y/L  ,  which  is  determined  by  two  "internal"  parameters, 
the  primary  peak  frequency  (which  gives  L0 ,  the  primary 
wave  length)  and  the  variance  which  are  routinely 
recorded  in  the  experiments.  Extension  of  this  model  to 
the  finite  depth  water  is  discussed.  The  controlling 
parameter  becomes  Ursell  number,  Ur’*H/k‘  h*  ,  which 
includes  the  depth  as  expected.  This  model  compares 
reasonably  well  with  the  experiment  data  in  the  energy 
containing  frequency  band.  At  higher  frequencies,  the 
deviation  becomes  generally  large.  An  alternative  of 
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using  the  stream  function  theory  to  simulate  the  spectrum 
shews  far  better  agreement  for  the  full  frequency  range. 
The  stream  function  theory  also  shows  good  agreement  with 
the  horizontal  and  vertical  velocity  spectra. 

4.2  A  Brief  Summary  of  the  Wallops  Spectral  Model 

Detailed  exposition  of  the  Wallops  spectral  model 
can  be  found  in  Huang  et  al.  (1981).  For  the  sake  of 
completeness,  a  short  description  will  be  given  as 
follows. 


The  form  of  the  Wallops  spectrum  is  based  on  a 
narrowband  Fourier  expansion  of  the  surface  wave  form 

l*Iah  cos  X„  (4-1) 

n 

with  X  is  the  phase  function,  and  an  is  the  amplitude  of 
each  Fourier  component  ne  simplest  form  of  this 
expansion  for  deep  water  waves  is  the  Stokes  wave, 

1  -  a  cos  X  +  1/2  a*  k,  cos  2X  (4-2) 

where  a  is  the  amplitude  of  the  foundamental  wave  and  k, 
is  the  wave  number.  Based  on  extensive  observations  in 
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the  laboratory,  Huang  et  al.  found  that  the  slope  of  the 
high  frequency  side  of  the  spectrum,  m,  can  be  determined 
by  the  relative  magnitudes  and  positions  of  the  primary 
and  the  sencond  harmonic  components.  Consequently, 

m  =  log(a,Va*)/log(2  <K/T. )  =  2  log  (  J2  rr§) /log  2  (4-3) 


With  the  high  frequency  part  determined,  the  whole 
spectrum  for  deep  water  waves  can  be  expressed  in  a 
similar  way  as  the  generalization  of  Pierson  and 
Moskowitz  (1964)  applied  to  Phillips'  saturation  range 
spectrum, 


e<pj'r 


(4-4) 


where  Se  denotes  the  wave  surface  soectrum  and 


r~-l 

von 


\ 


(4-5) 


in  which  ( )  is  the  gemma  function.  The  Wallops  spectrum 
approaches  Pierson-Moskowitz  spectrum  at  saturation 
stage  or  m  -*•  5.  It  also  eliminates  the  need  of  empirical 
coefficients,  all  the  parameters  are  internally 
determined,  and  the  spectrum  correctly  represents  the 


total  energy  content  through  the  definition  of  /S.  More 
detailed  discussion  of  this  model  can  be  found  in  the 
original  paper  (Huang  et  al.,  1981). 

4.3  Wallops  Spectrum  in  Finite  Depth  Water 

waves  in  water  of  finite  depth  behave  quite 
differently  from  the  deep  water  ones.  As  the  depth  of 
the  water  becomes  shallower,  the  effect  df  the  bottom 
will  be  felt  by  the  waves.  This  will  produce  important 
changes  in  the  wave  profile,  phase  velocity  and 
dispersion  relationship.  This  section  discusses  the 
modification  of  the  Wallops  spectrum  using  the  finite 
depth  Stokes  theory,  cnoidal  theory,  solitary  theory  and 
stream  function  theory.  While  the  former  three  theories 
predict  reasonable  agreement  near  the  main  peaks  for  a 
limited  kh  range,  the  stream  function  theory  shows  much 
better  agreement  over  the  whole  frequency  range  and  the 
range  of  relative  water  depth  of  this  experiment  (kh*0.2 
to  1.0).  The  derivation  of  each  case  follows. 


Finite  Depth  Stokes  Waves 
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A  profile  of  the  Stokes  wave  in  intermediate 
water  can  be  found  in  Stokes  (1880),  Bowden  (1958),  De 
(1955)  and  Wehausen  and  Laitone  (1960).  All  the  profiles 
are  essentially  equivalent.  Let's  use  the  expression  by 
Stokes  and  Bowden, 


b 

i 


a  ccs  X  + 


ai  k 


coth  kh  (!▼ 


1  OfVlk’  bh 


)cos  2x  (4-6) 


Then  following  the  same  approach  as  outlined  in  the 
previous  section,  we  can  write  the  slope  of  the  high 
frequency  range  of  the  spectrum  as 

m  »  2  log[  2TrJ  coth  kh  (1+  3  ^  /  log  2 

and 

a  -  _ ! _  (4-8 ) 

r  *«-•  \  1  *  0 1 

4.  -ar  P(  -3T  ) 

Comparing  the  expressions  of  Eqs.  (4-7)  and  (4-8)  with 
(4-3)  and  (4-5)  the  difference  is  obvious.  In  water  of 
finite  depth,  we  need  two  parameters  and  kh,  to  define 
m  and  (3,  The  values  of  m  and  <3  as  functions  of  §  and  kh 
are  shown  in  Figures  4.1  and  4.2.  As  kh  increases,  the 
influence  of  depth  decreases.  By  the  time  kh  reaches  3, 


ig.  4.1  Variation  of  m  and  5  based  on  Stokes 
finite  depth  wave  theory 
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both  m  and  essentially  assume  their  deep  water  values 
as  given  in  Huang  et  al.  (1981). 

Note  that  as  wave  travels  from  deep  water  to 

.  i 

. ifcermediate  depth,  evolution  occurs.  The  dotted  lines 
on  Fig.  4.1  shows  the  change  of  §  based  on  the  energy 
flux  conservation, 

c^E  *  constant  (4-9) 


and 


-  § /Itanh^kh  (1  +  2kh/sinh  2kh)]^  (4-10) 

where  denotes  the  deep  water  value  of  ^  ,  and  c,  is  the 
group  velocity. 

According  to  Stokes  expansion,  the  slope  can 
become  0  or  even  negative,  indicating  that  the  amplitude 
of  the  second  harmonic  is  higher  than  the  primary,  and 
the  Stokes  expansion  becomes  invalid.  In  fact,  long 
before  m  reaches  zero,  the  approximation  of  the  Stokes 
expansion  would  deteriorate  beyond  acceptable  limits. 
The  gradual  breakdown  of  the  Stokes  expansion  in  shallow 
water  is  also  reflected  in  the  values  shown  in  Fig. 


4.2.  Since  ft  is  derived  from  the  condition  of  correct 
total  energy  content,  a  finite  integrated  value  of  the 
whole  spectrum  is  necessary.  As  the  m  value  decreases 
below  1,  however,  the  Gamma-function  in  the  ft  expression 
has  logrithmic  singularities.  Consequently,  the  lower 
bound  for  the  applicability  of  the  Wallops  spectrum  will 
have  to  be  m>l.  An  alternative  approach  to  avoid  the 
singularities  is  to  define  a  multi-slope  spectral  model 
with  a  matching  frequency  at  a  few  times  of  the  primary 
frequency.  The  determination  of  this  matching  frequency 
need  more  comparison  of  experimental  data.  In  any  case, 
as  kh  decreases  to  about  unity,  we  need  to  use  other  wave 
models  for  modification. 

Cnoidai  and  Solitary  Waves 

As  the  waves  propagate  into  shallower  water 
(kh<l)  the  Stokes  expansions  eventually  break  down.  This 
can  be  seen  in  the  growth  of  the  harmonic  component 
relative  to  the  primary  waves  as  shown  by  Flick  et  al. 
(1981).  In  shallow  water,  the  cnoidai  wave  offers  a  much 
more  satisfactory  representation  of  the  waves.  The 
cnoidai  wave  model  was  developed  by  Korteweg  and  DeVries 
(1895)..  At  the  limits,  the  cnoidai  wave  approaches  the 
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solitary  wave  for  kh<<l,  and  sinusoidal  wave  for  kh>>l, 
which  is  not  exactly  (but  very  close  to)  the  Stokes  wave. 

The  surface  elevation  of  the  water  according  to 
the  cnoidal  wave  is 

7-  H  +  cn3  [9(A)  ,W})  (4-11) 

where  H  is  the  trough- to-crest  height,  Im'.n  is  the  depth 
of  the  water  measured  to  the  trough,  and  cn(&(-u),//J  is 
the  Jacobian  elliptic  function  with  the  phase  function, 

e-UO  »  [K(«  ) AOk(x-ct)  (4-12) 

and  the  cnoidal  parameter  ju , 

^[K (*)]a  «  (3Trl/4)  (H/kahJ )  (4-13) 

The  water  depth  to  the  trough,  is 

n  . 

--  (1  -  EU)/K (M))/M  -  1  (4-14) 

where  E(-w)  and  K(^<)  are  the  complete  elliptic  integrals 
of  the  first  and  second  kind  respectively. 
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For  the  limiting  case  of  x-~l,  K^)"**5  ,  then 

(-w),"]  =  sech'3 [  ( 3H/4h3  )yi  (x-ct)  J  (4-15)  « 

and  Iwin-^h,  so  that  the  surface  profile  becomes  I 

'l  *  h  +  H  seen2  C  ( 3H/ 4h'3  <x-ct))  (4-16)  I 

which  is  the  solitary  wave  solution.  The  limiting 

processes  adopted  here  is  very  critical  near  the  point  I 

For  example,  if  >-'  =  0.9959,  K ' M)»  6.0.  Therefore, 
for  practical  application,  a  wave  number  can  be  I 

introduced  to  make  the  wave  profile  of  finite  wave  | 

length,  Eq.  (4-15)  becomes 

0  =  h  +  H  seen2 [ ( 3H/4kJ  hJ )^k ( x-ct ) ]  (4-17) 

Both  cnoidal  and  solitary  wave  profiles  are  controlled  by  j 

a  single  parameter,  the  Ursell  number,  i.e., 

1 

Ur'  =  H/(k2hJ) 

The  Fourier  expansions  of  the  profiles  are 


a 

B 

a 
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a  X/H  «  ( - - — )  n  q’V  ( 1  -  q*n)  (4-19) 

3Ur ' 

for  the  cnoidal  wave  as  given  by  Cayley  (1895);  and 

a n/H  *  ( 1 — —  )/sinh[n7T/(3Ur  '  )^3  (4-20) 

3Ur ' 

for  the  solitary  waves  as  given  by  Gradshteyn  and  Ryzhik 
(1980).  The  quantity  g  in  Eg.  (4-19)  is  defined  as 

g  *  exp[  -II'k  ( l-xi)  /k  (M )  ]  (4-21) 

With  these  equations,  we  can  proceed  to  develop  the 
Wallops  spectral  model  in  shallow  water  using  either 
cnoidal  or  solitary  wave  as  the  basic  components  that 
make  up  the  wave  field. 

For  solitary  wave  case,  the  slope  of  the  high 
frequency  range  of  the  spectrum  is 

ms»  2  log  (a , /a  ,.)/log  2*2  log  cosh  M/  log  2  (4-22) 

with 


M  *  nr/  ( 3Ur '  ),/!l 


(4-23) 
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For  the  cnoidal  wave,  the  spectral  slope  is 


mc=  2  iog[  ( l-rq1)/2q]  /  log  2 


(4-24) 


For  both  these  cases,  the  {3  coef  f  ic ient  is  given  by  the 
general  expression 


(SJlS 


tr>  -t 


(4-25) 


where  c  is  the  phase  velocity  and  m  is  the  spectral  slope 
for  each  case.  For  most  studies,  the  phase  velocity  of 
the  Stokes  wave, 


c2=  (g/k)tanh  kh  (4-26) 


will  give  a  highly  accurate  answer  as  pointed  out  by  Bona 
et  al.  (1981).  with  this  expression,  the  parameter 
given  in  Eq.  (4-5)  could  be  used  for  both  the  cnoidal 
and  the  solitary  waves  provided  the  proper  m  values  as 
given  in  (4-22)  and  (4-24)  were  used. 

For  practical  applications,  it  is  more  convenient 
to  use  the  RMS  variance  as  the  wave  height  scale 
(especially  for  field  data).  The  new  Ursell  number 
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becomes 


Ur  *  T/Z/  (k2h3)  =  Ur'  »)*% 


(4-27) 


CjVi. 

where  7  is  the  RMS  elevation  with  respect  to  the  mean 
water  level.  For  a  sinusoidal  wave  train,  the  mean  is 
zero,  hence 


H/<2/2) 


(4-28) 


For  the  solitary  waves,  the  mean  water  level  is 


ff  =  h  + 


H 


^j^tanh  [ir(3UrV4)  ] 


( 4-25) 


and 


[(  n  -  T  )a  ]  '/a 


-[***  («-*wa-  ,4-30) 


where  Q  *  TT( 3Ur ' /4 )v*,  and  n'*mean  elevation. 


The  value  of  (4-30)  is  given  in  Fig.  4.3.  The  special 
significance  of  this  quantity  is  that  it  represents  the 


r  LO 
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ratio  of  Ur'/Ur,  Eq.  (4-27).  Since  most  data  are 
expressed  in  terms  of  Ur,  Eq.  (4-30)  provided  the  link 
between  observations  and  the  model  results.  The 
corresponding  values  for  the  cnoidal  waves  are 

H  {1  -  [l-EUO/K U)]/*}  (4-31) 

and 

—  I/n 

h“  =  I  ~(2-4-u)n'l-  r*\ 

(4-32) 

The  value  of  Eq.  (4-32)  is  plotted  on  Fig.  4.3  also. 

Stream  Function  Wave  Theory 

Stream  function  wave  theory  developed  by  Dean 
(1965)  was  based  on  nolinear  least  square  fitting  of  the 
governing  equation  and  the  boundary  conditions  of  a  wave 
field.  The  solution  is  in  the  form  of 

MM 

V(x,z)-zl/T  +2L  X(n)sinh[^(h+z)]cos(-2n?x) 

(4-33) 
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and 

*)  =T  ^j/L  -T/L  2L  X(n)sinh[^p;'(h+,1)  Icosfny^x) 

(4-34) 

The  solution  of  Eqs .  (4-33)  and  (4-34)  satisfy  the 

governing  Laplace  equation  and  bottom  and  free  surface 
kinematic  conditions.  L  and  X(n)  are  numerically 
calculated  to  best  fit  the  dynamic  free  surface 
condition.  The  detailed  treatment  is  presented  in  Dean 
(1965,  1974).  using  the  procedure,  by  specifying  the 

wave  period,  water  depth  and  wave  height,  the  surface 
elevation  can  be  obtained.  Applying  Fourier 
decomposition,  the  spectrum  of  each  wave  can  be  obtained, 
which  will  be  termed  "stream  function  spectrum".  The 
corresponding  Wallops  m  and  can  be  calculated  from  the 
first  and  second  harmonic  values  and  the  variance  of  the 
elevation  record.  For  practical  applications, 
computation  on  a  high  speed  digital  computer  is 
necessary. 


4.4  Velocity  Spectra 
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For  engineering  applications  such  as  wave  force 
analysis  for  marine  structures,  sediment  transport 
mechanisms  and  nearshore  circulation,  the  ability  of 
correctly  predict  the  velocity  field  is  important.  The 
velocity  field  under  a  nonlinear  wave  field  is  a  complex 
one.  This  section  only  presents  the  horizontal  and 
vertical  velocity  spectra  measured  in  the  sediment 
experiment.  Comparison  of  the  measurement  with  the 
linear  transformation  from  wave  spectra  and  the  stream 
function  velocity  spectra  are  discussed. 

Based  on  the  linear  wave  theory,  the  velocity 
spectral  transformation  can  be  expressed  as, 

Su('J‘)*tTccsh  k (h+z)/sinh  khTSeCT)  (4-35) 


and 


Sw(<T)-[<3sinh  k(h+z)/sinh  kh]aSe(3‘)  (4-36) 

where  Suf  Sw  and  Se  are  the  spectra  of  the  horizontal 
velocity,  the  vertical  velocity  and  the  wave  surface 
elevation,  respectively.  On  the  other  hand,  the 


velocities  based  on  stream  function  wave  theory  can  be 
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written  as 

u  =  -  21  X(n)  (27Tn/L)cosh£  ZJL  n(h+z)  ]cos  (—px)  (4-37) 

n=.|  <—  *— 


and 


w  *  21  X(n)  (2im/L)sinh[-^-n(h+z)  JsinC^px)  (438) 

The  u  and  w  profiles  can  then  be  processed  by  Fourier 
decomposition  to  create  the  velocity  spectra.  The 
comparison  with  experimental  data  is  shown  in  next 
section . 

4.5  Comparison  with  Observations 

Having  constructed  the  spectral  models  for 
intermediate  and  shallow  water  waves,  they  are  tested 
against  the  observations.  As  a  first  step,  the  values  of 
the  spectral  slope  in  the  high  frequency  range  is  plotted 
as  a  function  of  Ur  in  Fig.  4.4.  Since  the  Stokes  wave 
depends  on  %  and  kh  in  a  mutually  unrelated  way,  the 
slope  value  should  be  represented  by  a  family  of  curves 
each  for  a  specific  §  value.  For  the  present  experiment 
data,  §  varies  from  0.005  to  0.025,  the  corresponding 
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Stokes  curves  are  shown.  The  numerical  results. of  the 
stream  function  calculation  also  show  nonunique  Vocation 
on  Ur.  Since  wave  length  is  a  variable  to  be  determined 
in  the  stream  function  calculation,  the  family  of  curves 
are  plotted  with  different  H/h  ratios.  The  numerical 
procedure  also  applied  to  each  individual  cases  of  the 
sediment  experiment  (Chapter  3)  and  plotted  with  symbol 
"F".  Wave  data  from  the  drift  experiment  (Chapter  2), 
Flick  et  al.'s  shoaling  experiment  (1981)  and  Thornton's 
field  experiment  (1979)  were  also  analyzed.  Thornton's 
data  is  in  spectral  form;  the  spectral  slopes  near  the 
energy  containing  part  are  used.  It  is  seen  that  for 
Ur>0.3,  cnoidal  and  solitary  results  are  nearly  identical 
and  show  reasonable  agreement  with  the  data  for  Ur>l.  As 
Ur  decreases  (deeper  water),  they  started  deviating  from 
the  experimental  results.  In  this  intermediate  depth 
range,  Stokes  result  shows  better  agreement.  Only  stream 
function  shows  good  agreement  for  the  whole  Ur  range 
tested. 


The  full  spectra  of  surface  elevation  for  the 
sediment  experiment  are  shown  in  the  top  plots  of  Fig. 
4.5.  For  clarity,  only  the  Wallops  spectra  based  on 
cnoidal  theory  and  stream  function  curves  are  plotted 
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Fig.  4.5  Comparison  of  1,  u  and  w  spectra  with 
various  theories 


Fig.  4.5  (continued) 
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Fig.  4.5  (continued) 
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Pig.  4.5  (continued) 
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along  with  the  experimental  results.  In  these  plots, 
only  the  spectral  peak  values  are  shown.  The  agreement 
over  the  higher  energy  components  is  generally  good  for 
the  Wallops  model,  but  for  the  full  range  of  frequency, 
direct  calculation  of  the  stream  function  wave  form  shows 
far  better  agreement.  It  is  noted  that  the  spectra  based 
on  the  cnoidal  and  solitary  theories  will  predict  a  much 
steeper  spectrum  beyond  the  second  or  the  third  harmonic 
as  implied  in  Fig.  4.6.  The  figure  compares  the 
theoretical  amplitude  harmonics  with  those  of  the 
measured  wave  profiles.  It  is  obvious  that  after  the 
second  harmonic,  the  theories  underestimate  the  harmonics 
and  therfore,  the  spectral  components.  The  middle  plots 
in  Fig.  4.5  are  the  horizontal  velocity  spectra  at  the 
elevation  of  z/h=-0.6  to  -0.7.  Better  representation  of 
the  stream  function  calculation  over  the  linear 
transformation  is  obvious.  The  bottom  plots  are  the 
vertical  velocity  spectra  at  the  same  depth.  Stream 
function  calculation,  again,  shows  very  good  agreement. 

Fig.  4.7  plots  the  variance  ratios  of  the  wave 
and  velocity  spectra  based  on  the  stream  function 
calculations  and  the  linear  transformation.  For  the  wave 
spectra  (plot  a),  the  mean  value  of  the  ratios  of  the 
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Fig.  4.6  Harmonic  amplitudes  Vs.  kh 
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stream  function  calculated  variances  to  those  of  the 
measurements  is  0.95  with  standard  deviation  of  0.37. 
The  variance  of  the  Wallops  model  for  m>l  are  also 
plotted  on  the  same  figure.  The  ratio  mean  and  deviation 
are  0.92  and  0.46,  respectively.  The  Wallops  spectral 
model  seems  to  represent  the  total  energy  adquately. 
Similar  plots  for  the  horizontal  and  the  vertical 
velocity  specta  are  given  in  Fig.  4.7  b  and  c, 
respectively. 

For  the  horizontal  velocity  spectral  variance 
(plot  b),  linear  transformation  performed  reasonably  well 
in  deeper  water  (Ur<l)  but  has  a  tendency  of 
over-prediction  as  water  becomes  shallower.  The  stream 
function  calculation  remains  good  and  stable  for  the 
range  tested.  There  are  some  cases  with  Ursell  number 
between  1.3  to  2.0  show  large  deviation.  Although  the 
definite  reasons  are  uncertain,  one  possible  explanation 
is  the  reflection  effects.  If  the  measuring  location 
happens  to  be  at  a  quasi-antinodal  point,  the  horizontal 
velocity  will  decrease  due  to  the  standing  wave 
contribution  at  the  antinode  (thus  a  higher  theory  to 
measurement  ratio).  This  is  partially  revealed  by  the  T 
and  w  variance  ratios,  most  of  which  are.  less  than  unity 


in  that  Ur  range.  More  extensive  studies  are  necessary 
to  clarify  this  point.  The  ratio  mean  and  deviation  are 
1.12  and  0.51  for  stream  function  theory,  and  1.53  and 
0.78  for  linear  transformation. 

Comparison  of  the  vertical  velocity  variance 
ratios  (plot  c)  is  similar  to  the  horizontal  one  except 
that  the  vertical  ratios  are  more  scattered.  Since  the 
magnitude  of  the  vertical  velocity  is  smaller  than  that 
of  the  horizontal  velocity,  any  factors  contributed  to 
the  deviation  such  as  turbulence  or  the  instrument  noise 
will  show  a  more  prominent  effect  on  the  vertical 
velocity.  The  ratio  mean  and  deviation  are  0.84  and  0.48 
for  the  stream  function  calculations,  and  1.91  and  1.20 
for  linear  transformation. 

In  summary,  the  Wallops  spectral  model  can  be 
modified  for  finite  depth  application.  The  model  will 
give  accurate  energy  content  and  good  agreement  in  the 
vicinity  of  the  primary  peaks  where  energy  are 
concentrated.  The  applicability  of  Stokes'  modification 
is  about  kh>0.75.  For  shallower  water  cases,  the  cnoiaal 
or  solitary  theory  is  more  appropriate.  For  m<l,  the 
Wallops  model  predicted  an  infinite  energy  content. 
Truncated  spectrum  or  multi-slope  model  need  to  be 
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considered. 

An  alternative  to  the  spectral  model  is  to 
calculate  the  spectra  from  the  wave  profiles  generated  by 
the  stream  function  wave  theory.  From  the  test  cases 
presented,  it  is  found  that  the  wave  and  velocity  spectra 
of  the  stream  function  generated  waves  show  excellent 
agreement  with  the  measurements. 


CHAPTER  5  CONCLUSION  AND  RECOMMENDATIONS 

Flow  properties  including  both  mean  and 
instantaneous  velocities  as  well  as  sediment  suspension 
processes  at  the  breaking  zone  were  studied  in  the 
laboratory.  Analysis  were  performed  in  the 
characteristics  of  the  mean  drift  velocity  and  the 
spectral  presentations  of  the  fluctuating  components  of 
water  surface  elevation  and  flow  velocities.  The 
mechanics  of  sediment  suspension  under  the  influence  of 
wave  motion  is  critically  examined. 

Drift  velocity  profile  at  breaking  point  shows  an 
onshore  motion  near  the  surface  and  close  to  the  bottom. 
In  the  middle  colume,  the  flow  is  always  offshore.  The 
magnitude  of  the  current  is  more  uniform  than  those  occur 
under  nonbreaking  waves.  The  shapes  of  the  profiles 
under  plunging  and  speilling  waves  are  similar. 
Comparison  of  the  vertically  averaged  return  flow 
velocity  with  four  drift  velocity  theories  was  conducted. 
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Reasonable  agreement  was  obtained. 

Nonlinear  effect  of  breaking  caused  harmonic 
generation  of  the  monochromatic  incident  waves.  The  wave 
form  was  analyzed  by  spectral  analysis.  Comparison  with 
the  nonlinear  stream  function  wave  theory  shows  good 
agreement  both  in  terms  of  the  individual  spectrum  and 
the  total  variance  calculation.  The  difference  is 
generally  under  15  percent.  The  Wallops  spectral  model 
modified  for  finite  water  application  was  proposed.  The 
model  gives  correct  variance  estimate  and  emphesizes  the 
simulation  of  the  higher  energy  containing  harmonics. 
The  instantaneous  velocity  measurement  was  compared  with 
the  stream  function  wave  theory,  the  difference  is 
generally  under  20  percent. 

The  sediment  suspension  can  be  represented  by  a 
turbulent  diffusion  model.  The  diffusivity  is  related  to 
the  wave  parameters,  with  wave  height  as  the  length  scale 
and  the  vertical  velocity  as  the  velocity  scale  for  high 
agitation  condition.  For  low  agitation,  the  local 
vertical  excursion  amplitude  is  the  proper  length  scale. 
Fall  velocity  reduction  in  an  oscillating  flow  field  and 
the  bottom  sediment  grading  effects  were  also  included  in 
the  model.  The  model  gives  good  agreement  both  in  the 
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concentration  profile  and  the  magnitude. 

For  future  studies  of  this  problem,  the  author 
suggests  the  following  directions, 

1.  Instrumentation:  The  instruments  for  sediment 

concentration  measurement  and  velocity  measurement  should 
be  improved.  The  frequency  response  should  be  in  the 
same  range  for  all  the  devices.  The  sensing  volume  of 
the  velocity  meter  need  to  be  kept  small  to  increase  the 
depth  resolution.  To  reduce  the  intermittency  of  the 
concentration  data,  larger  sensing  volume  to  enhance 
chance  of  particles  passing  through  is  necessary.  The 
tradeoff  of  the  signal  quality  and  the  frequency  response 
need  careful  studies. 

2.  Tracing  of  neutrally  buoyant  particles  to  obtain 
Lagrangian  drift  is  subjected  to  the  processor's 
judgement  and  introduces  error.  Other  procedure  to 
reduce  human  interpretation  is  desired  for  a  more 
accurate  measurement. 

4.  Different  concepts  on  suspension  mechanisms  that 
emphasize  the  vortex  trapping  has  been  offered  recently 
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(Nielsen,  1981;  Tooby  et  al,  1977).  Longuet-Higgins 
(1980)  offered  a  numerical  model  of  vertex  formation  and 
shedding.  It  is  always  useful  to  have  a  different 
perspective  to  one  problem.  The  result  may  enhance  our 
knowledge 

5.  Wave  theories  considering  the  slope  effects  had  not 
been  included  in  this  preliminary  study.  For  a  more 
refined  model,  the  slope  effects  need  to  be  included. 
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Appendix  Lagrangian  Mass  Transport  by  Stream 
Wave  Theory 

Dalrymple  (1976)  employed  a  Eulerian  reference 
frame  to  evaluate  wave  mass  transport  using  Dean's  stream 
function  wave  theory  (Dean,  1965,  1974).  He  also 
extended  the  theory  to  include  the  cases  with  the  linear 
shear  current  superimposed  on  the  wave  field  (Dalrymple, 
1973,  1974).  In  this  appendix,  a  method  of  evaluating 

Lagrangian  drift  velocity  by  following  particle  motion 
along  a  streamline  is  developed.  The  drift  velocity 
profiles  are  similar  to  those  of  Stokes  (1847).  The 
total  mass  transport  is  the  same  as  evaluated  by  Eulerian 
method,  which  serves  as  a  check  in  the  validity  of  the 
numerical  scheme. 

Theory 

The  mass  transport,  M  is  defined  as 


M»  J^u  ( z ,  t )  dzdt 


(1-1) 
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or  changing  the  order  of  integration 

M-f /_”  [  1/T  jfju(  z  '  , t )dt ]dz '  (1-2) 

where  u(z',t)  is  the  Lagrangian  velocity  following  a 
particle.  The  bracketed  term  in  Ea.  (1-2)  is  known  as 
the  mass  transport  vel'  ity  or  drift  velocity,  U.  Dean's 
stream  fuction  theory  states  that  the  stream  function  of 
a  wave  field  can  be  expressed  as 

M'fx  ,z )  =zL/T  ■*J£x(n)sinh[*fc?(h+z)  ’cosf-^x)  (1-3) 

L.  L-— 

where  coordinate  system  moves  with  wave  celerity  to 
render  steadiness  of  the  wave  motion.  The  particle 
velocity  in  the  wave  field  can  be  written  as 

u(x,z)«-l-  7-7  XfnJcoshf—^ (h-rz )  ]cos(^5^x)  (1-4) 

The  elevation  z'  of  a  particular  streamline,  'V',  can  be 
obtained  by  iteration  from  the  following  expression 

z'(x)*T»V/L  -  T/L  3L  X(n ) sinh[ ~p(h+z '  )  ]cos ( ~-x )  ( 1-5 ) 

Using  Eq.  (1-4)  to  evaluate  M,  a  transformation  of 
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variable  from  t  to  x  is  necessary, 

Ms^Clou{t'z'  )«3tdz'-^jf(Lu(x,z'  )|~~dxdz’  (1-6) 

The  Jacobian  in  Eq.  (1-6)  is  J=2t/3x.  In  stream 

function  theory,  the  coordinate  system  is  moving  with 
wave  celerity,  C.  That  is,  at  any  wave  phase 
© ■x-Ct«constant.  Therefore,  J«l/C.  Eq.  (1-6)  becomes 

fLu(x,z')/C  dxdz'  (1-7) 

Using  this  equation,  however,  we  have  to  calculate  u 
above  water  surface  for  the  upper  layers.  This  is 
physically  wrong  especially  for  nonlinear  waves  which 
have  long  troughes  below  mean  water  level.  In  the 
following  paragraph,  a  different  approach  to  evaluate  M 
is  presented.  The  mass  transport  velocity  is  calculated 
by  following  a  streamline,  therefore,  particles  always 
stay  in  the  fluid  region. 

By  using  Sq.  (1-5)  to  obtain  the  elevation,  z'. 


of  a  streamline,  particle  velocity  along  the  streamline 
can  be  calculated  by 


u(x,z'  )=  2.  X(n  )cosh[-^p[h+z '  )  3cos(i£px)  (1-8) 

The  wave  length  is  then  divided  into  N  segments  of  ax 
each.  The  small  time  interval,  at  that  the  particle  will 
travel  a  distance  ax  can  be  computed.  The  net  particle 
displacement  as  well  as  the  net  velocity  is  obtained  if 
the  computation  is  carried  out  over  one  wave  cycie,  i.e., 

U(z')= - Z.  u(x.-  ,Z’-)4t;  (1-9) 

ZAt;  <-  1  1 

)  1 

where  ati^^x^/u-,  z'  denotes  the  mean  elevation  of  the 
streamline  followed.  M  is  then  calculated  by  streamline 
followed.  M  is  then  calculated  by 

M-f  i  lOfT'j  )-U(^,)](z;-z^,)  (1-10) 

Mass  Transport 

To  check  the  validity  of  the  above  numerical 
scheme,  the  mass  transport  computation  using  Eq.  (1-10) 
is  compared  with  Dalrymple’s  Eulerian  computation. 
Column  6  and  7  of  Table  1.1  list  the  mass  transport 
calculated  by  these  two  different  methods  (col.  6  - 

present}  col.  7  -  Oalrymple).  The  orders  of  wave  theory 
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Table  1.1 

Lagrangian  vs.  Eulerian  Mass  Transport 
by  Stream  Function  Theory 


1 

2 

3 

4 

5 

6 

7 

CASE 

H 

T 

h 

ala-  k 

kh 

M/f  Bu, 

(m) 

(s) 

(m) 

(m/s) 

(m/s) 

(m/s) 

A-l 

10.24 

10 

30.48 

1.492 

1.246 

16.63 

16.54 

c-i 

5.32 

10 

30.48 

0.203 

1.369 

2.47 

2.47 

C-2 

6.10 

10 

30.48 

0.261 

1.363 

3.22 

3.21 

C-3 

9.10 

10 

30.48 

0.574 

1.332 

7.01 

7.00 

C-4 

10.63 

10 

30.43 

0.765 

1.313 

9.25 

9.23 

C-5 

16.14 

10 

30.43 

1.656 

1.233 

17.78 

17.63 

C-6 

21.29 

10 

30.48 

- 

- 

- 

28.59 

B-l 

1.92 

10 

3.05 

0.058 

0.306 

0.414 

0.412 

D-l 

0.60 

10 

3.05 

0.006 

0.346 

0.071 

0.070 

D-2 

0.91 

10 

3.05 

0.014 

0.336 

0.145 

0.144 

D-3 

1.91 

10 

3.05 

0.024 

0.327 

0.222 

0.221 

D-4 

1.52 

10 

3.05 

0.038 

0.317 

0.315 

0.313 

D-S 

1.77 

10 

3.05 

0.050 

0.310 

0.381 

0.378 

D-6 

2.38 

10 

3.05 

0.087 

0.297 

0.478 

0.460 

calculation  were  17  pnd  13  for  h/gT=0. 00311  (cases  A  and 
C),  and  h/gT*0.G311  (cases  B  and  D) ,  respectively.  A 
wave  length  is  divided  into  72  equally  spaced  segments. 
As  shown  in  the  table,  only  negligible  difference 
present . 

The  profiles  of  Lagrangian  mass  transport  for  a 
wide  range  of  kh  values  are  given  in  Fig.  1.1.  The 
profiles  of  Stokes  theory  are  also  given  for  comparison. 
For  deepwater  cases,  both  profiles  resemble  each  other. 
When  kh  decrease,  the  difference  oetween  these  two 
profiles  becomes  larger.  The  Stokes  prof ile ( becomes  more 
and  more  uniform  over  the  depth.  The  stream  function 
profile,  on  the  other  hand,  shows  a  very  strong  forward 
drift  near  the  free  surface.  Since  both  theories  do  not 
consider  the  boundary  effects,  no-slip  condition  at  the 
bottom  boundary  is  not  satisfied.  Table  1.2  summarize 
the  conditions  and  the  results  presented  in  Fig.  1.1. 
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Table  1.2 

Comparison  between  Mass  Transport  Computed  by  Stokes  Theory 
and  Stream  Function  Theory 


1 

2 

3 

4 

5 

6 

7 

CASE 

H 

T 

h 

aV  k 

kh 

M 

M 

(m) 

(s) 

<m) 

(m/s) 

a 

1«rk(fh) 

a*7k  (fh) 

A-l 

15.24 

10 

30.48 

1.492 

1.246 

0.390 

0.474 

3-1 

6.29 

10 

3.05 

0.058 

0.306 

2.411 

5.505 

E-l 

0.050 

0.44 

0.15 

0.155 

2.658 

0.155 

0.190 

E-2 

0.094 

0.70 

0.15 

0.157 

1.217 

0.368 

0.490 

E-3 

0.111 

0.99 

0.15 

0.098 

0.764 

0.613 

1.017 

E-4 

0.116 

1.40 

0.15 

0.050 

0.503 

1.000 

2.140 

E-5 

0.118 

2.21 

0.15 

0.020 

0.299 

1.953 

5.758 

E-6 

0.118 

3.12 

0.15 

0.009 

0.205 

3.258 

12.064 

E-7 

0.113 

4.42 

0.15 

0.005 

0.156 

12.000 

20.712 

Column  6:  Stream  function  theory 
Column  7:  Stokes  theory 
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